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ABSTRACT 

We argue that most strong intervening metal absorption line systems, where the rest equivalent width of the 
Mg IIA2796 line is > 0.5A, are interstellar material in, and outflowing from, star-forming disks. We show 
that a version of the Kennicutt-Schmidt law is readily obtained if the Mg II equivalent widths are interpreted as 
kinematic broadening from absorbing gas in outflowing winds originating from star-forming galaxies. Taking 
a phenomenological approach and using a set of observational constraints available for star-forming galaxies, 
we are able to account for the density distribution of strong Mg II absorbers over cosmic time. The association 
of intervening material with star-forming disks naturally explains the metallicity and dust content of strong 
Mg II systems as well as their high H I column densities, and does not require the advection of metals from 
compact star-forming regions into the galaxy halos to account for the observations. We find that galaxies with 
a broad range of luminosities can give rise to absorption of a given rest-equivalent width, and discuss possible 
observational strategies to better quantify true galaxy-absorber associations and further test our model. We 
show that the redshift evolution in the density of absorbers closely tracks the star formation history of the 
universe and that strong intervening systems can be used to directly probe the physics of both bright and faint 
galaxies over a broad redshift range. In particular, in its simplest form, our model suggests that many of 
the statistical properties of star-forming galaxies and their associated outflows have not evolved significantly 
since z ~ 2. By identifying strong intervening systems with galaxy disks and quantifying a version of the 
Kennicutt-Schmidt law that applies to them, a new probe of the interstellar medium is found which provides 
complementary information to that obtained through emission studies of galaxies. Implications of our results 
for galaxy feedback and enrichment of the intergalactic medium are discussed. 
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1. INTRODUCTION 

Quasar spectra often reveal the presence of interven- 
ing metal absorption systems that span the entire observ- 
able redshift range b etween the observer and the source 
(Burbidge et al. 1966). A possible association of such ab- 
sorption systems w i th inte rvening galaxy disks was first sug- 
gested by Wagoner ( 1967). An alternative explanation associ- 
ating such systems with large (compared t o the size of visible 
disks) halos of galaxies was considered by Bahcall & Spitzer 
(1969), who paved the way to the use of intervening ab- 
sorption sys tems as probes of the large scale environment 
of galaxies (IWeymann et al.l Q978). Despite enormous ef- 
forts over the last three decades to determine whether ab- 
sorption arises in galaxy disks or halos, the question re- 
mains largely unresolved. This is particularly true for the 
Mg II AA2796,2803 absorption line systems; these absorbers 
can be detected at relatively low redshifts (z > 0.2) in the opti- 
cal spectra of quasars, so understanding their origin has been 
a more tractable problem for modern instrumentation, as we 
now discuss. 

Initial follow-up imaging and spectroscopy of galax- 
ies in the fields of quasars known to show Mg II ab- 
sorption confirmed the association of strong (with rest- 
equivalent width Wo > lA) intervening systems with galax- 
ies whose luminosities are in the range 
dBergeron|[l986t ISteidell [T995b . Nevertheless, a clear pic- 
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ture connecting absorbers and galax i es has yet to emerge 
dYorkl | 1982t iBergeron & Boisselfl991h iBechtold & Ellingsonl 
| 1992t Isteidel et alJI 1994 119971 12002 [Tripp & Bowenl l2005l: 
iKacprzak et alJl2007t IChen & Tinker! 120081) . Much progress 
has been made in characterizing the properties of absorp- 
tion line systems from studying the spectra of lar ge numbers 
of QSOs discovered in wide-area sky surveys dNestor et alj 
HOOI IProchter et all 120061 lYork et all 12009) . Following 
on from the earlie st analysis of the r e dshift distribution 
of absorption lines (iLanzetta et alJl!987t iTytler et all 119871: 
Sargent et alj!988l) it is now well documented that the cosmic 
density of absorbers per unit redshift and unit rest-equivalent 
width, as a funct ion of Wn, d 2 N I dzdWn, follows an exp onen- 
tial distribution dNestor et al. 2005; Pro chter et al.ll2006l) . but 
how this distribution is related to the association of absorption 
lines with galaxies, is unknown. Recent st atistical analyses 
reveal that strong M g II systems are dusty (I York et al J 120061: 
Menard et al. 2008) and that the dust-to-gas ratio of moder- 
ate redshift systems is consistent with local inte rstellar ma- 
terial (ISM) values (M enard & Chelouchdl2009L but may be 
substantially lower at higher redshifts dPetitjean et al.lll996l) ). 
This suggests that the gas composition of such systems is 
not considerably different fro m solar values (lYork et al.ll2006t 
Menard & Chel ouchel 12009b . There are recent indications 
that strong systems are associated with star formation activ- 
ity (iBouche et al.ll2006L 12007b IWild et"aTll2007t Izibetti et al.1 
l2007t IMenard et al. 2009) yet the nature of this connection is 
not well understood. 

Despite the recent progress in the phenomenology of in- 
tervening absorption systems, little is known with con- 
fidence about the physical state and origin of the gas 
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dChurchiU et alj|2005h IChelouche etaTI 120081 and references 
therein). Different models at t ribute t he absorption to either 
galact ic winds [e.g., Bouche (2008); Oppenhei mer & Pavel 
(2008)], or to cool clumps condensi ng out from the virial- 
ized g a seous halos of ga l axies l e.g.. IMo & Miralda-Escudel 
( 1996); lMaller & Bul lock (2004); Kauf mann et all (120081)1 re^ 
lated, perhaps, to infallin g material dTinker & Chenl 120081; 
iKeres & Hernquisll 120091) . Alternative explanations re- 
late some s ystems to gravitationally bound minihalos 
(iKeoner et al j(l999t ISternberg et al.ll2002t iGnat & Sternberg! 
2004) or to galaxy disks (Wagoner 1967; Bowen et al. 



1995; Prochaska & Wolfe 1997; Charlton & Churchill 1996; 
Prochaska & Wolfe! 1 1 998t IZwaan et al.l 120051 12008k While 



some models are qualitatively consistent with s o me ob - 
servational constraint s le.g.. ISrianand & Khard dl994l) ; 
IChelouche eTall (120081) : iTinker & Chenl (120081) 1. many of the 
recently discovered statistical properties of Mg II systems 
have not been addressed by theoretical works, and the fol- 
lowing questions remain: What is the origin of strong metal 
systems? What sets their observational signatures? How do 
they relate to galactic and intergalactic medium (IGM) pro- 
cesses? In this we attempt to shed light on some of those 
questions by establishing a firmer connection between strong 
intervening systems and star forming galactic disks. 

The paper is organized as follows: section 2 qualita- 
tively argues for a Kennicutt-Schmidt like star formation law 
( Schm idl[T959l: lKennicutdfl98l fl998l) for strong intervening 
Mg II systems. A more quantitative approach is taken in §3 
where a phenomenological model is formulated, which links 
the properties of star-forming galaxies to those of absorbers. 
The results of our model fits to the observations are outlined 
in §4. We discuss the implications of our results for the study 
of strong intervening systems and galaxy disks across cosmic 
time in §5. A summary follows in §6. 

2. A KENNICUTT-SCHMIDT LAW FOR STRONG MG II ABSORBERS 

Perhaps the best known property of star forming galaxy 
disks is that they obey the Kennicutt-Schmidt la w for star for- 
mation dSchmidtlll959l: lKennicutj[l989l 119981) whereby the 
star formation rate per unit area scales with the mean column 
density of the disk to some powefl The underlying physi- 
cal mechanism behind this phenomenological scaling is not 
fully understood and considerable theoretical and observa- 
tional effort has been invested in identifying its origin (|Silk| 
1997|:lKravts^l2003blKrumholz & McKed2005l:lBigiel et alj 
2008[ iLerov et all 120081; iKrumholz et all 120091: b enzel et alj 
2O10tlGnedin & Kravtsovll2010l) . 



Recent spectroscopic observations reveal that many star 
forming galaxi es give rise to outflows of cool ga s from 
their centers (iHeckman et al.l 119901 l2000t iMartinl [2005) 
Ruoke et alJ|2005h IVeilleux et al]l2005t iTremonti et al.N2007t: 
Weiner et alj|2009t : iRubin et alj|2009t) . Furthermore, a rela- 
tion exis ts between the outflow velocity and the star forma - 
tion rate (lMartinl2005l: I Weiner et al.ll2009t IRubin etal]l2009l) . 
The physical nature of such outflows and their prevalence 
in the galaxy population as a whole, are currently poorly 
constrained. Nevertheless, it has been suggested that such 
outflows can provide the necessary negative feedback ef- 
fects that determine and mainta in the observed star forma- 
tion rates observed in galaxies dStrickland & Stevens! [T999I : 
LShaviv & DekeMX B: Murr ay et allm)¥ . and as a means for 

3 We note, however, possible deviatio ns from a simple powerlaw at low 
column densities and star formation rates iBigiel et al. 2008); see also §5.8 
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FIG. 1 . — A schematic view of a wind launched from a galactic disk (here 
we neglect inclination effects). The length of the arrows does not reflect 
the extent of the outflow above the disk and is for illustrative purposes only. 
The cross-section of the wind may be smaller than that of the H I disk but 
is unlikely to be much smaller than that of the UV emitting disk (see text). 
The spectrum of background luminous objects (whose light is not signifi- 
cantly attenuated by extinction), such as quasars and y-ray bursts (GRBs) 
will show absorption systems whose velocity dispersion is determined by the 
wind kinematics on the approaching and the receding side of the disk. UV 
observations of galaxies in which star forming regions are used as the back- 
ground UV source probe only the approaching side of the disk outflow. We 
neglect contributions to the absorption signal from infalling material. Incli- 
nation effects (not shown ) may be i mportant and are incorporated into our 
model using data from Martin (2005). 



enrich ing the IGM dDalcantonl 120071 : lOppenhei mer & Pavel 
2008, and refer ences therein). 

As mentioned above, there is mounting evidence that strong 
Mg II absorbers (Wo > lA) are connected to star formation 
activity and, as such, co uld be related to the gaseous out- 
flows associated with it (iNorman et al.l Il996t iProchter et alJ 
120061; iBouche et al]|2007l) . Here we argue that strong Mg II 
absorbers are, in fact, star forming galactic disks and their as- 
sociated outflows. To this end we note two observationally 
robust scalings: the relation between the (geometric) mean 
H I column density of strong M g II absorbers and Wo such 
that dMenard & Chel ouche 20Ql), 



JVhi = (3 ±0.5) x 10 iy W c 



19,,,1.7±0.3 







cm 



(1) 



and the relation between the (maximum) velocity of galactic 
winds, v w , and the star formation rate (SFR, denoted by T) 
such that 



1M yr 



-l 



0.35±0.05 



km s 



(2) 



Here, t, w = 1 for Na I lines dMartinl2005l) and t, w < 3 for Mg II 
lines, as statistically estimated from the stacked spectra of a 
sample of star-forming DEEP2 galaxies dWeiner et al.ll2009b . 
To clarify we note that % w > 1 does not necessarily mean that 
Mg II ions travel faster than Na I ions and that it may be that 
observational effects, su ch as optical depth or partial filling by 
emission are important dMartin & Bouc he 2009). 

Here we assume, with little justification at this stage, that 
outflow kinematics are responsible for the line broadening of 
strong Mg II systems. In particular, we consider a case in 
which the rest equivalent width of strong Mg II troughs and 
the velocities of outflows are related by 



W ~2 



(^21 ) 
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This proportionality between the kinematic broadening and 
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Wo is justified for strong Mg II systems whose absorption 
troughs are saturated, and Wo mirrors the velocity dispersion 
of the gas. Since v w is a measure of the outflow velocity to- 
wards the center of a galaxy, we add the proportionality factor 
of 2 to account for a background object sightline intercepting 
both the receding and advancing sides of a symmetric outflow 
(see Fig. [TJ. This model neglects the contribution o f infalling 
material fe-g., cool filam ents or galactic fountains (Bregman 
1980; Brooks et al. 2009)] whose covering factor is expected 
to be small dSato et al.l 12009) and the velocity spread mod- 
erate. It also neglects contribution to the velocity field from 
galactic rotation which is expect ed to be small for t hin disks, 
as determined from current data (Zwaan et al. 2008). Inclina- 
tion effects will be treated in §3. 

Combining the above, seemingly unrelated, relations for 
iVffl(Wo) an d using equation[5] we obtain an expression 
relating the SFR to the H I column density, with the following 
form: 



T/L (lCr 10 M o yr- 1 L~ r . 
10~ 1 10° 




T~(0.1±0.03): 



HI 



V 10 20 c 



,-2 



1.7±0.3 



M yr 
(4) 

This relationship between the SFR of a Mg II host galaxy 
and the H I column density of the absorbing gas, qualitatively 
has the same functional form as the Kennicut t -Schmidt law 
for star formation in gala xies (Kennicutil ll998t iBouche et ail 
I2007t iGenzel et al.l l2010). An analo gous scaling is obtained 
if, instead of using equation Q] one uses th e relation be- 
tween the dust column density and Wo found bv lMenard et al.l 
d2008lfl 

Let us now crudely estimate a Kennicutt-Schmidt rela- 
tion from equation |4] we recall that Wo > l A systems ar e 
thought to be associated with L* galaxies (Steidel 1995), 
whose disk radius we take to be R e f[ ~ l Okpc [qualita- 
tively consistent with the relation reported by ICourteau et al.l 
(2007)]. An order-of- magnitude estimate for the star for- 
mation surface density is therefore £* = Y/(tcR^ s ) ~ 3 x 
1 O- 4 <i 7 2 ± o°' 3 (^/3)- 2 - 9±O ' 4 M yr- 1 kpc" 2 (where N m , 20 = 
Nfti/10 2Q cm~ 2 ). We hence obtain a surp risingly good ma tch 
with the galactic Kennicutt-Schmidt law (Kennicuttl ll998l) . It 
is important to note that the uncertainties on the proportional- 
ity factor in this relation are large: for example, if the typical 
galactic disk radius giving rise to strong systems is only 3 kpc 
large (as would be the case if the hosts of most strong sys- 
tems are fainter; see below) then the Kennicutt-Schmidt law 
obeyed by absorbers implies a higher star formation rate for 
the same H I column density compared to galaxies. We return 
to these important issues when considering a more quantita- 
tive version of this relation in §5.8. 

Obtaining the form of the Kennicutt-Schmidt law by com- 
bining pure absorber and galaxy phenomenology raises the 
interesting possibility that (a) sight-lines through most strong 
Mg II absorbers prob e star forming disks fas oppo sed to e.g., 
condensed halo gas; (Mo & Miralda-Escude 1996)] and that 
(b) the broadening of Mg II troughs is largely due to disk wind 
kinematics rather than, for example, the velocity dispersion 
of absorbing clouds in galactic halos or the rotation speed of 
gas in galactic disks. This intimate absorber-disk association 

4 IMenard ef"aT] <2008[) find that the reddening E(B - V) due to intervening 



FIG. 2. — The observationally derived cumulative probabilities [C(x) = 
Jn P(.x')dx J ] used in this paper. C (Wp|Y) (solid lines for several SFR values) 
was obtained from Martin 12005) (after conversion to Wo; see text). Clearly, 
an order of magnitude range in Wo may be obtained for any given SFR bin 
-2 9±0 4 [ see also[Martin (2005) where this holds also for narrow bin sizes]. Gener- 

_ i ally, however, smaller SFR result in lower outflow velocities, hence Wo val- 
ues. C(T/L) (dashed lin e) shows the specific SFR distribution as determined 
from James et al. 1 2008). Clearly, a broad SFR range can be obtained for any 
given galaxy. 

also naturally explains the large H I column densities asso- 
ciated with strong intervening systems, which are comparable 
to those observed through disk s. It also explains the high d ust- 
to-gas ratios in those systems (Men ard & Chelouchel2009l) . as 
well as their association with star forming regions. 

3. THE MODEL 

If most strong Mg II systems probe galaxy disks then their 
statistical properties should be accounted for by those of 
galaxies. In particular, the exponential density distribution 
of such absorption systems as a function of Wo, d 2 N/dzdWo 
should be describable by the following expressiorjj 



systems' dust and Wo are related such that E(B — V) °= W ' 



Assuming 



d 2 N{W ,z) 
dzdW 



(1 



#o E(z) Jl, 



dL<j> z (L)o z (L)P z (W \L) 



the reddening is a proxy for the dust column and a constant dust-to-gas ratio, 
a completely analogous relation to equation [TJis obtained. 



(5) 

where E(z) — \J£Im{^ +z) 3 +^a, <!> z (L) is the redshift- 
dependent luminosity function, L m ; n (L m ax) is the minimum 
(maximum) luminosity of galaxies that give rise to strong 
Mg II systems, and O z (L) is the luminosity and redshift- 
dependent geometrical cross-section for detecting a Mg II ab- 
sorption around a galaxy of luminosity L. P(Wo|L) if the prob- 
ability of detecting a system with a rest equivalent width Wo 
around a galaxy of luminosity Z0. This probability may be 
written as 

P Z (W Q \L) = J dYP z (W Q \Y)P z (Y\L), (6) 

where P Z (Y\L) is the probability of having a certain SFR, T, 
for galaxy of luminosity L, and P,(Wo|Y) is the probability 
of having a certain Wo - i.e., outflow velocity - for a galaxy 
with a star formation rate T. For simplicity we assume that 
there are no hidden correlations among the various variables 
apart from those explicitly stated here. Whether this holds in 
realistic systems is unclear, but assessing the degree to which 

5 We assume a flat universe with £1a = 0.7, Q.m = 0.3, and Hq = 
70km s- 1 Mpc" 1 . 

6 Unless otherwise specified, the following normalization holds: 
fP(x\y)dx=l 
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any hidden correlations may change our results is beyond the 
scope of this paper. 

The model defined above seeks to connect the observable 
properties of absorbers to those of their host galaxies. Un- 
fortunately, the properties of galaxies and absorption line sys- 
tems are well measured only over particular redshift intervals, 
and these intervals do not always overlap. The properties of 
galaxies are best understood at very low redshift, but the char- 
acteristics of UV absorption line systems have been deter- 
mined largely at intermediate and high redshift. For example, 
as we shall see below, faint galaxies have a potentially impor- 
tant contribution to the cosmic density of strong systems and 
their properties are known only locally. We can therefore pro- 
ceed in one of two ways: we can (a) extrapolate the local prop- 
erties of (faint) galaxies to higher redshifts, or (b) extrapolate 
the properties of absorbers to low redshifts. The first approach 
requires the extrapolation of the locally determined galaxy 
luminosity function, O z =q(L), P z=0 (Y\L), and P z =o(Wo\Y) to 
higher z. In contrast, the second approach requires only that 
we extrapolate d 2 N(Wo)/dzdWo to low redshifts. Further- 
more, with the installation of the Cosmic Origin Spectrograph 
(COS) on the Hubble Space Telescope (HST), the properties 
of low redshift absorption lines are likely to be determined (at 
least for W) > lA systems) in the near future. Conversely, 
many of the properties of faint galaxies at high redshifts are 
likely to remain uncertain for some time. 

With the above in mind, we restrict our analysis to z ~ 
and will henceforth drop the z— dependence from our derived 
values and distributions (we return to the effects of redshift 
evolution in §5.5). Specifically, we take a Schechter form for 
the luminosity function of galaxies that give rise to metal ab- 
sorption systems 

m dL= fL r^y e -^f (?) 

We adopt the parameterization of lJames et al.l (120081) who car- 
ried out a census of star-forming galaxies in the local uni- 
verse, and found: a ~ — 1.4, 0* ~ 3 x 10~ 3 Mpc~ 3 (L* ~ 
1.4 x 1O 1O L0, with Lr.j being the solar luminosity). We note 
that this particular parameterization for the luminosity func- 
tion of galaxies is one out of many, and somewhat differ- 
ent, parameter values are reported in the literature depending 
on the selection criteria and filte r used (Marzke et al. Il998t 
iZwaan et alJl200U iNorberg et afll200a iBlanton et al.H2003l) . 
As we do not know whether all star-forming disk galaxies give 
rise to Mg II systems, we define fi to be the fraction of galax- 
ies that cause Mg II absorption and allow it to be smaller than 
unity. We neglect a possible dependence of fi on the galaxy 
luminosity. We take L max — > °° implying that the model con- 
siders all galaxies with luminosity, L > L m ; n as contributing 
to the population of Mg II absorbers. It turns out that, due to 
the exponentially small number of very luminous galaxies, the 
particular L max chosen has little effect on the predicted num- 
ber density of Wo > lA absorbers. That said, the particular 
choice of L max may affect model predictions concerning the 
number density of the strongest systems, which are rare and 
their statistics rather poorly determined. 

The cross-section for Mg II absorption from galaxy disks, 
d(L) = 7tRMgii{L) 2 [with /?Mgii(i) being the effective radius 
of the absorbing region], is unknown, as is its dependence on 
galaxy luminosity. We assume that this cross-section scales 
with that of the H I disk, Ohi, as deduced from 21 cm emission 



down to a column density limit of > 10 19 cm 2 ( IZwaan et alj 
120051) . We therefore define 

o(L) = &aa(L*)(j£) r (8) 

with 7 being of order unity, which is qualitatively consis- 
tent with the obse rved local scaling of H I disk cross-sections 
[see Fig. 10 of IZwaan et al.l (120051) . which shows a con- 
stant probability for detecting absorbers through disks span- 
ning a r ange of luminosities] . For infrared galaxy disks 
y~ 0.7 dCourteauet all 120071) while in the optical y ~ 0.6 
dde Jon g & Lacev 2000). Interestingly, a similar scaling 
has been recently reported between the effective radius of 
Mg II absorbing e nvelopes of galaxies and their luminosity 
dChen et al.ll2010l) . In equation [8] ^ a stands for the relative 
size of the Mg II strong absorption cross-section compared to 
the H I disk (assumed to be luminosity independent). This 
factor is currently poorly constrained and will be determined 
below from a fit of our model predictions to available absorp- 
tion line data. 

We next define the proba bilities P(Y\L ) and P (Wo | Y) based 
on data from the literature. iJames et"alj d2008l see their Fig. 
3) found no clear trend between the specific SFR (i.e., SFR 
per unit luminosity), as determined from Ha luminosities, 
and galaxy luminosity; we determine P(Y/L) from their data, 
and present its cumulative version in Figure [2] The proba- 
bility of having a certain rest equivalent width for the Mg II 
line per give n star f ormation rate, P{Wq\Y), was determined 
from iMartinl ((2005) who measured the outflow velocity in 
Na I absorption lines as a function of the SFR. In doing so we 
assume that outflow kinematics, as observed along a sight- 
line to a backgro und object, traces global outflow proper- 
ties ((Martin 2005), and neglect possible complications related 
to t he effects of partial cove ring of the wind over the disk 
area lMartin & Bouchel (12009). The translation from the out- 
flow velocity of Na I lines to Wo, using equation [3j requires 
knowledge of £L (see §2) whose value € [ 1,10] ( Martirj |2b05h 
lTremonti etaiTl2007t iWeiner et al] 120091 see also §2) but is 
left here as a free parameter and will be determined from a fit 
of equation [5] to the available absorption line data. We hav e 
used data pertaining to the high SFR bins of Martin (2005), 
resulting in the best characterization of P(Wo = § w v w |Y) for 
high SFRs (see her Fig. 6 for T > 100 M yr" 1 ). We then as- 
sume that the deduced P(Wq\Y) also holds for lower SFR up 
to a scaling factor such that the maximum velocity, hence Wo, 
in each star formation bin is °< Y? with j8 ~ 0.3 dMartinl2005l: 
Wein er et al.l 120091 see our Fig. This is consistent with 
the IMartinl d2005l) results and is, presently, the only viable as- 
sumption given the quality of the data. To summarize, our 
model thus far relies on observable properties of star-forming 
galaxies with only two quantities, £, w and £, a , both poorly con- 
strained by observations and left as free parameters whose 
values will be determined by a fit to absorption line data. 

It is, currently, not well determined how inclination affects 
the observable properties of a galaxy outflow and any attempt 
to directly account for it in our model is subject to large un- 
certainties. While it is unlikely that an outflow is spherically 
symmetric, the details are obscure both observationally and 
theoretically [for a very r ecent work provid ing new insights 
into inclination issues see IChen etalJd2Ol0h l. That said, our 
model is unlikely to be considerably affected by such uncer- 
tainties since we use the observationally determined statisti- 
cal properties of outflows from spectroscopic surveys of star 
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FIG . 3. — The d 2 N(z _ 0) / dzdWo distribution with the normalization taken 
from Prochter et al. 1 2006) and the uncertainty estimated to be a factor of 
~ 3 for Wo ~ 1, and somewhat larger for larger Wo values (see §4). The 
shape of the d 2 N(z = 0)/dzdW curve was determined from Ne stor et"aTI 
12005) by extrapolating their high redshift data (see their Fig. 10) to lower 
redshifts using their non-evolutionary curve models. Good agreement be- 
tween the model and the data is shown in black dashed line and its parame- 
ters are: S = a + y= -1.4 + 0.7 = -0.7, j3 =0.3, = 1.6, K Mg n(i*) = 
4Okpc,0* = 3 x 10~ 3 Mpc~ 3 , L rain = 10~ 2 L* (model A; see table 1). Also 
shown are a few additional models, all of which are consistent with the data 
given the uncertainties, in colored curves (unless otherwise specified param- 
eter values are identical to model A): £ w = 2 (green), !; w = 1.1 (orange), 
£ w = 1.1, 8 = -0.4 (magneta), P = 0.25 (cyan), R Mg n(i*) = 50kpc (blue), 
and T] = 0.1, § w = 3.6 (red; see §4.1). Clearly, the origin of all strong systems 
with Wo > 0.7A at z — is consistent with local star forming disks. Diamonds 
mark our evolutionary model predictions for the signal at z. = (see §5.5). 



forming galaxies. Such studies probe galaxies with a broad 
range of inclinations to our sight line, hence the effects of in- 
clination already appears in the data and enter our model, by 
definition, via P(W \T). 

We predict the number density of absorption line systems 
with a given equivalent width, d 2 N(Wo,z = 0)/dzdWo from 
equation[5]using a Monte Carlo scheme: we first pick a galaxy 
whose luminosity, L, follows from the luminosity function. 
Given L, we use P(T\L) to statistically determine its star for- 
mation rate. The following step statistically determines Wo 
from P(Wo|T). At this point, we have fully specified the 
statistical relation between galaxy and absorbing properties. 
This scheme is iterated to create a mock catalog of galax- 
ies and their absorbing characteristics from which the statisti- 
cal properties of the predicted absorbers may be deduced and 
compared to observations. The size of the sample is deter- 
mined by the requirement that good statistics are realized for 
the strongest Mg II absorbers, which are inherently rare. 

3.1. The Kinematic Properties of Outflows Across Disks 

The properties of Mg II winds across galaxy disks are 
poorly determined and only a characteristic velocity, v w , over 
the entire UV emitting disk surface is measured. Whether or 
not outflows are thermally expanding or radiatively driven, it 
is likely that their kinematics will depend on their proximity 
to luminous or energetic sources in the disk such as star form- 
ing re gions, dense stellar environme nts, or active galactic nu- 
clei (Martin 2006; Fuiita et al. 2009, and references therein). 
These are not distributed uniformly over t he galaxy but rather 
cluster in the inner ann uli of the disk dLerov et al.l 120081; 
Rahman & Murrav 2010, and references therein). To allow 
for the (likely) possibility that gas kinematics varies across 
the disk surface of galaxies, we consider a general velocity 



profile of the form: 

Vw(p) = &(p)v» ■= + &" v » (9) 

with p being the impact parameter to the host's center and 
e < is a powerlaw index (with E, w having the usual defi- 
nition). We define the (projected) core radius within which 
the velocity profile flattens as some fraction of the radius of 
the Mg II absorbing disk such that, Pq(L) = T]RMgii(L) with 
77 assumed constant. We choose e = —0.5 so that on scales 
p 3> po, an inverse square-root law is recovered mimicking, 
for example, the r adial dependence of the e scape velocity 
from a point mass (IChelouche & N etzer 2005) or the termi- 
nal velocity of a radiation pre ssure driven cloud from a point 
source of a given luminosity (Chelouche & Netzer 200JJ). On 
scales p < po galactic structure becomes important and the 
approximation of a point source breaks down leading to a flat- 
tening of the velocity profile. Having very little information 
about the dynamics of galactic outflows, this is the simplest 
form for the velocity profile that is consistent with the data 
(see §4. 1 for further details The inclusion of this part of the 
model in our Monte Carlo calculations is straightforward and 
requires us to randomly pick the impact parameter p in addi- 
tion to the host luminosity in order to predict Wo for a specific 
sightline. We present results for this extension to the model in 
§4.1. 

At this point, our model now relies on three free parame- 
ters: £,a(L*) [or RMgii(L*)], and r\ all of which are poorly 
determined by observations. In what follows we shall restrict 
our analysis to models in which 77 is either 1 (a constant 
wind velocity over the disk surface) or tj = 0. 1 which roughly 
corresponds to the star-forming region of galaxies being en- 
closed within their inner few kpc. 

4. RESULTS 

When comparing our model predictions to the data it is 
important to note that many of the parameters in equation 
[5] (e.g., n z , ov) are constrained by observations of galax- 
ies at low redshifts and that their high-z analogs may, in 
principle, have different values. Instead of attempting to 
predict how galaxy properties may have changed with red- 
shift, we instead extrapolate d 2 N/dzdWo to redshifts of zero. 
As the statistics of strong Mg II systems are only loosely 
determined at present times, we have used higher redshift 
data to estimate d 2 N/dzdWo at low-z. In p articular, we 
used th e non-evolutionary curves in Figure 10 of lNestor et al. 
(2005) to estimate the relative normalization, i.e., the shape 
of d 2 N/dzdWo, at different rest equivalent width slices used 
by the authors over the range 0.6 < Wo < 3. 5 A. The overall 
normalization of d 2 N/dzdWo at z — was chosen such that 
J^dW d 2 N(W ,z = 0)/dzdW Q = dN{W > lk,z = 0)/dz = 

0.05 

^ao3' which is consistent with the extrapolation s to z — 
based on the best-fit models of iProchter et al.l (120061 see their 
Table 4) to their high redshift data set. The quoted un- 
certainty qualitatively brackets the range in values predicted 

7 Not all radial forms are allowed: cuspy velocity profiles for which the 
velocity of the wind diverges as p — > [as in the case, v w (p) °= p e for which 
rj <@C 1] would be in disagreement with the data since, in this case, high Wo 
values in the d 2 N / dzdWo distribution would be determined by low luminos- 
ity objects rather than high luminosity ones. In this case, model predictions 
will take on a powerlaw form at large Wo values reflecting the geometrical 
cross-section for fast outflows rather than the exponential decay of the lumi- 
nosity function (see below). 
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by the different extrapolation methods used in iProchter et alj 
( 120061) and the measured value of O.I6+QQ5 reported by 
IChurchilll d200U see his Fig. 7) for z ~ 0.05 systems. The 
number statistics of stronger systems (Wo > 2A) is less se- 
cure at low redshifts, and our uncertainty for Wo > 2A bins 
was chosen to include the factor between the predictions of 
the non-evolutiona ry models and the evolutionary curves of 
iNestor et alj d2005L compare their Figs. 10 and 13) for z = 0. 
We neglect sample incompleteness effects arising from dust 
extinction at the high Wo end of the distribution but note that 
these can result in an under-estimation of dN(Wn > 3A,z — 
1 ) /dz by up to - 20% dMenard et al.ll2008l) . 

The numerical evaluation of equation [5] for the case of a 
constant wind velocity over the disk surface, and its compar- 
ison to the observed number density of Mg II lines, as ex- 
trapolated to z = 0, is shown in Figure [3] (The solution for 
a model in which the wind varies as a function of p, with 
77 < 1, is discussed in §4.1.) In this case, the fitting scheme 
involves a search in a two parameter space that is spanned by 
% w and £, a that are rather poorly determined observationally 
[the former parameter enters implicitly through the definition 
of P(W)|I") (Fig. |2} while the latter enters explicitly in equa- 
tion |5). Here we have taken a = — 1.4, y = 0.7 (we define 
8 = a +J as this is the only relevant parameter given the form 
of Eq. [5] hence 8 = —0.7), and j3 = 0.3. While somewhat 
uncertain, these values are relatively well constrained obser- 
vationally and the sensitivity of the model to the small uncer- 
tainties in those parameters will be further investigated below. 
A good fit to the (extrapolated) d 2 N(z — 0)/<?z<?W) data is 
obtained for t, w = 1.6 and for an effective cross-section for 
Mg II absorption around an L* galaxy whose projected radius 
Ru g n{L*) = ^o{L*)/k = ^Z, a o m {L*)/7l satisfies: 



R U ga(L*) =s 40 



1 dN(W > lA) 



,1/2 



0.05 



dz 



0* 



1-1/2" 



fl 



-1/2 



3 x 10- J Mpc- 



(10) 



kpc. 



IZwaan et alj (120051 see their Fig. 18) find that the typical (pro- 
jected) H I radius of an L* galaxy at z — 0, as detected down 
to a column density limit of 3 x 10 19 cm~ 2 is ~ 35 kpc, i.e., 
of order 7?Mgii for /z, = 1, indicating that % a ~ 1. Such scales 
are also consistent with the size of the Mil ky Way disk mea- 
sured at similarly low H I column densities (Kalbe rla & Pedes! 
120081) . 

We examined how galaxies of different luminosities con- 
tribute to the estimated values of d 2 N / dzdW n. We show 
the results i n Figure ffl If , as po inted out by Steid el et al.1 
(IT991 fl997l) : iKeeney et al.1 (120061) . galaxies with L <C 0.1L*, 
are seen to give rise to strong absorption, then absorbers with 
Wo > 0.5A can indeed be fully accounted for by star-forming 
galaxy disks. Whether L <C 10~ 2 L* galaxies give rise to 
strong metal absorption systems is unclear, since our model 
does not require it. Furthermore, at such low luminosities, 
galaxy disks become H I poor a nd do not seem to follow the 
relation assumed by equation[H](Zwaan et al. 2005|). 

It is important to note that our model does not reproduce 
the number density of weak absorbers with Wo <C 0.5A. As 
can be seen in Figure [4] at these equivalent widths, there is a 
clear discrepancy between the number density of Mg II lines 
predicted for all galaxy luminosities in our model, and the 
observed number density of Mg II systems. This difference 
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FIG. 4. — The contribution of galaxies of different luminosities to the oc- 
currence of Mg II intervening systems. The contribution of all galaxy disks 
with luminosity L > L m j„ to the d 2 N(z, = 0)/dzdWo signal is shown for sev- 
er al values of L^ m d enoted next to each dashed curve. The best-fit model 
of INesto r et al. 12005) to their data (solid line) is scaled up by 20% to al- 
low an easy comparison with the model. Clearly, weak absorbers with 
Wo < 0.5A cannot be accounted for by the model even when the contribu- 
tion from very faint galaxies, with L -C 10~ 2 L*, is taken into account. The 
residual signal, after the subtraction of the L mm = 10~ 3 L* model from the 
data, may be interpreted as the weak absorbers' signal contributing to the to- 
tal d N(z = 0)/dzdWo distribution, and is not accounted for by our model 
(dotted line). The signal for strong absorbers is dominated by relatively lu- 
minous star forming, galaxies. 

suggests that the physical origins of wea k Mg II systems ma y 
be different from that proposed herein (Richter et al. 2009). 
This conclusion is related to and further strengthened by the 
fact that Wo ~ 0.5A is also the valu e at which a break is seen 
in the d 2 N(z > 0)/dzdW Q function (INestor et al.ll2005l) . 

Consider now the sensitivity of the predictions for d 2 N(z > 
Q)/dzdWo to the various parameters of the model. In partic- 
ular, £, a (or /?Mgii ) is essentially unknown, as is £, w . [One 
may obtai n indirect statistical information regarding E, w from 
the data of Weiner et al. (2009) yet its interpretation in terms 
of ^ w is not unique; see §5.4]. We find that our model pre- 
dictions change in the following way: varying E, w changes the 
shape of the predicted curve such that larger \ w results in a 
more shallow decline at the high- Wo end. This is due to fainter 
and more numerous galaxies contributing to the signal at the 
high-Wo end. Varying ^ a [or /?Mgii(i*)] merely changes the 
normalization of the models but maintains their shape unal- 
tered. The same behavior applies for /z,. In addition to the 
above parameters, we have also considered the effect of mea- 
surement uncertainties on the observationally deduced values 
used by our model and held fixed in the above analysis. For 
example, adopting smaller values of j3 result in more low lu- 
minosity objects giving rise to stronger lines giving rise to a 
flatter d 2 N(z > 0)/dzdWo distribution. Likewise, increasing 
8 results in a smaller contribution to the number density from 
low luminosity objects (since the product of their density and 
cross-section - i.e., their opacity - is smaller), creating a flatter 
d 2 N(z > 0)/ dzdWo curve for weak lines. 

In order to quantify the degeneracies inherent in the model, 
and to constrain the physics most relevant to forming strong 
intervening systems in galaxy disks, we have considered a 
range of values for all the relevant parameters in our model: 
20 < R Mg n(L*) < 60kpc (0.3 < \ a < 2.3), -0.7 < 8 < 
-0.1, 0.25 < j3 < 0.4, 10~ 3 <L min /L* < 0.1, and 1 < % w < 3. 
In this multi-parameter space, we have calculated a grid of 
models and evaluated their agreement with the data using % 2 
statistics. Examples for a few such models, and their agree- 
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FIG. 5. — The rejection confidence contours for our models. The multi-parameter space has been divided to subplots, each of which refers to two parameters 
among the five and the results for the remaining parameters have been marginalized over. Each point gives the mean rejection probability of the marginalized 
model, as determined from reduced x 2 statistics, x 2 - The mean value of x 2 at eacn point (not shown) is the geometric mean of the reduced x 2 of all models being 
marginalized over. The 60% contour level (see labels) correspond to a Xv — 1- B should be emphasized, however, that even in regions whose (mean) rejection 
probability is relatively high, there may be individual models providing a very good fit to the data. Clearly, model predictions are in good agreement with the 
data for the various parameter ranges considered here. Model predictions are particularly sensitive to the values of R-Mgll{L*) and t; w with the other parameters 
having a lesser effect on the quality of the fit. 



ment with the data are shown in Figure [3] (dotted curves). 
Presenting the regions in phase space for which the fits are 
acceptable, we show in Figure [5] various contour plots as a 
function of two parameters at a time. The quantity plotted is 
the mean significance with which models may be rejected. It 
should be emphasized that, in each inset, we have marginal- 
ized over the remaining parameters of the model such that 
the geometric mean of the % 2 values has been evaluated at 
any point and the mean significance accordingly calculated. 
We also note that, due to our marginalization, it is possible, 
in principle, to find very good agreement between specific 
models and the data also in parameter space regions where 
the mean significance is low. Clearly, the models are in good 
agreement with the data for the range of parameters discussed 
above. For example, models for which the size of the Mg II 
disk of an L* galaxy, is greater than ~ 30kpc (for fi = l and 
with the disk size of other galaxies following from eq. |8j, can 
explain the observed d 2 N/dzdWo distribution, given the un- 
certainties. To conclude, the fit is sensitive to £, a and <jj w , and 
to a lesser extent on the values of 8, /3, and L m ; n within the 
observationally motivated range of values. This also means 
that our model is, effectively sensitive only to % w , and tj 
and therefore provides a < 3 parameter fit to the data. 

It is evident from Fig. [5] that good fits are obtained for 
% w > I. In particular, the model wherein the outflow veloc- 
ity field across the disk is uniform, suggests that 1 < £, w < 2. 
Therefore, our model predicts that the Mg II transitions ex- 
tend to somewhat higher velocities than the Na I lines. This 



is reasonable given the higher magnesium abundance with re- 
spect to sodium (while neglecting ionizing corrections), and 
the stronger oscillator strengths of the lines in the former case. 
Alternatively, the apparently higher velocities for the Mg II 
tra nsitions could be a purely observational effect, as discussed 
in dMartin & Bouche 2009). Our conclusion that % w > 1 is 
further supported by the difference in t he Mg II line profile s 
seen towards the centers of galaxies by Weiner et al. (2009), 
compared to the Na I profiles observed by Martin (2005|). We 
return to this point in §5.4. 

If most star-forming galaxies brighter than 10~ 2 L* give rise 
to strong Mg II absorption then Figure [5] suggests that 35 < 
^Mgii(i*) < 60kpc for the studied paramet er space. Such 
sizes are comparable to those deduced by iKacprzak et al. 
(2008, see their eq. 3) given the occurrence of strong sys- 
tems and based on a very different set of arguments. Our 
deduced sizes are al so consistent with the recent findings of 
Mena rd et al.l ((2009) who show that the O II luminosity func- 
tion, which traces star-forming regions in the universe, can 
be fully accounted for by the number density of strong Mg II 
systems that surround star-forming galaxies on similar scales. 
Interestingly, such sizes have also been found to characterize 
the H I disks of L* galaxies, as measured for the Westerbork 
HI Survey of Irregular an d Spiral Galaxies (WH ISP) s ample 
( Ivan der Hulst et al.ll200lh . In particular, IZwaan et all (120051 
see their Fig. 18) found a projected H I radius of L* galaxies 
to be in the range 25-50 kpc. This implies that, in our model, 
~ 1. This conclusion is consistent with the findings of 
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iBowen et al.l dl995l see their Fig. 29) who showed that the 
projected radius of the H I disk for an L* galaxy is > 25 kpc, 
and that the size of the Mg II absorption cross-section is of 
at least a comparable size. This means that the cross-section 
for strong Mg II absorption is much larger than that of the 
optical and UV emitting disks. This makes outflows visible 
for studies in which the galaxy itself acts as a background 
source against which Mg II systems that are associated with 
the ga laxy are detectable (Tremonti et al. 2007; Wein er et alj 
2009). 

The match between H I disk sizes a nd the pred i cted R mvU is 
also consistent with properties of the iRao et alJ (2006) sam- 
ple of absorbers for which the typical H I column densi- 
ties that are associated with Wo > 0.5A systems are of order 
> 10 19 cm~ 2 . Such H I column densities are somewhat be- 
low the limiting column densities of galaxies in the WHISP 
sample, which is ~ 3 x 10 19 cm~ 2 dZwaan et al. 2005). Given 
that the apparent disk sizes increase with increasing sensi- 
tivity of the observations (so that lower gas columns can be 
probed), and extrapolating the results of IZwaan et al.l (2005, 
see their Fig. 18), we find that the projected radii of L* galax- 
ies down to a column density of 10 19 cm~ 2 is > 40 kpc, in line 
with our findings. We note that about 50%, of Mg II systems 
with 0.5 < Wo < lA have H I column densities smaller than 
5 x 10 18 cm dRao et al.l 120061) i.e., considerably below the 
detection limit of the WHISP sample. This implies that the 
cross-section for dilute low column density gas, as appropri- 
ate for W) > 0.5A systems, is essentially unknown and may be 
larger than that which is currently de duced from low sensitiv- 
ity 21 cm emission studies [see e.g.. lBraun & Thilkeri J200l 
who find evidence for an extended envelope of low column 
density gas on scales > 50 kpc around M31]. 

Further evidence that strong Mg II systems are probing 
disks comes from the fact that the coherence size of W) ~ 0.5 
absorbers, as estimated from observations of lensed quasars, 
is found to be 3^ kpc (Ellis on et al.l 120 04"). In our model, 
such systems predominantly arise in low luminosity galaxies 
with 10~ 3 L* < L < 10~ 2 L* (see Fig. 0]and below). As such, 
our model implies a coherence length of order the disk radius, 
i.e., R Mgll (L/L*y/ 2 ~ 3.5 - 8kpc (for/ L = 1), a value that is 
consistent with observations. 

4. 1 . The Kinematic Properties of Outflows Across Disks 

In §3.1 we defined the projected core radius of the disk, 
within which the outflow velocity approaches a constant 
value, to be po(L) = rjRMgii(L). In our model for the kine- 
matics of outflows across disks (Eq. [9]), we adopt 77 = 0.1 
for no other reason than the fact that, in this case, po coin- 
cides with the optical, star-forming size of the disk. Running 
Monte Carlo simulations and including the effect of random 
sight lines through the disk surface, we find that it is possi- 
ble to fit the data well if £, w — 3.6. This means that, for an 
L* galaxy, the core in the velocity profile occurs for Po(L*) = 
T]RMgii(L*) ~ 4 kpc, while the maximum velocity of the Mg II 
lines is > 3 times higher than that measured for the Na I 
lines. For galaxies with T ~ 1OM yr _1 , this model predicts 
maximum wind velocities of > 250 x 10 03 > 500km s" 1 to- 
ward their cores. Note, however, that the mean velocity (and 
conversely Wo), as measured for a statistical sample of such 
galaxies, will be lower since only a fraction of all objects 
would launch outflows from their centers that travel at the 
maximum velocity according to P(Wq\T) (see also §5.4). 
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FIG. 6. — Predictions for d 2 N '/ dzdWo for several values of and rj: mod- 
els with ~ Tj -0 - 5 and with all other para meters held fixed are in good 
agreement with the data of Nestor et al. (2005, dashed red curve) in the range 
0.5 < Wo < 4A. Nevertheless, for larger Wo, predictions are sensitive to the 
value of v w . In particular, larger v w , and correspondingly smaller rj, result 
in a powerlaw decrease rather than an exponential decay in the number of 
large Wo systems. This results from the number statistics of such systems 
being dominated by relatively faint star-forming galaxies for which quasar 
sightlines pass close to their centers rather than the exponential tail of the 
galaxy luminosity function. Note the sharp rise of the observed d N /dzdWo 
for Wo < 0.5A resulting from the unique population of weak systems, which 
is not accounted for by the model. 

The above solution is, however, non-unique, and equally 
good fits may be obtained, for example, when taking % w = 6 
and rj = 0.03 (while holding all other parameters fixed). In 
particular, the current data quality for d 2 N /dzdWo, as extrap- 
olated to z ~ 0, and for large values of Wo, provides little con- 
straints on the allowed rj ranges, and that ^ w ~ 77 ~ 05 
traces a good solution curve in phase space with all other 
parameters held fixed and given the current data; see figure 
[6] However, different models have different predictions for 
the number statistics of the strongest systems. In particu- 
lar, in models for which v w is higher (and rj correspondingly 
smaller), the abundance of the strongest systems is dominated 
by sightlines passing close to the centers of relatively faint 
galaxies (see above). Currently, the d 2 N/dzdWo data set 
alone cannot be used to reliably constrain physical models for 
galactic outflows. As we discuss in §5.4, however, there are 
additional data that may be used to shed light on the kinemat- 
ics of galactic outflows. 

5. DISCUSSION 

A Kennicutt-Schmidt law for strong absorbers is consistent 
with the notion that most strong Mg II systems seen in the 
spectra of background quasars arise in disks of intervening 
star-forming galaxies. We have shown that the density distri- 
bution of strong Mg II absorbers, d 2 N(z = 0)/ dzdWo, can be 
fully accounted for if galaxy disks with the same properties 
as those observed in low-redshift galaxies are responsible for 
the absorption. Our model implies that (a) the Mg II absorp- 
tion cross-section for a galaxy is comparable to that of the H I 
cross-section of nearby galaxy disks, and that (b) the veloc- 
ity of the outflowing gas that gives rise to intervening Mg II 
absorption is consistent, up to a factor of order unity, with 
the outflow kinematics measured for Na I and Mg II absorp- 
tion troughs towards the centers of local star-forming galax- 
ies. The notion that strong Mg II systems are material in and 
outflowing fro m galaxy disk s is als o consistent with the re- 
cent findings of iBernet et al.l d2008l) indicating that the mag- 
netic fields associated with strong systems are comparable to 
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galactic values. 

The association of strong Mg II systems with galaxy disks 
is also physically appealing as it naturally resolves several 
outstanding theoretical problems in the field: for example, 
if strong Mg II systems are interstellar material in extended 
disks then the problem of cloud formation, con finement, and 
stabil ity in dilute intergalactic environments dBinnev et alj 
2009) is, clearly, alleviated. Furthermore, there is no need to 
transport cool metal-rich material from galaxy centers to large 
scales while avoiding over-ionization and mixing. In fact, in 
our model, the physics of strong intervening Mg II systems is 
exactly that of the ISM, which is (at least initially) gravita- 
tionally bound and enriched by stars in its environment. 

Below, we discuss a few implications of our results. 

5.1. Absorber-Galaxy Associations 

As demonstrated above, a broad range of galaxy luminosi- 
ties gives rise to strong Mg II absorption. For Wo > 0.5A sys- 
tems, the galaxies which give rise to the absorption have lu- 
minosities > 10~ 2 L*. For example, the model shown in Fig- 
ure |4] predicts that there is a 50% chance that the hosts of 
Wo ~ 0.7A Mg II systems are fainter than 0. 1L*. It may there- 
fore be challenging to identify the galaxy which directly hosts 
the absorber, since it could be outshone by the nearby quasar. 
In such cases only bright hosts or neighboring galaxies that 
are clustered around the true host may be detected (provided 
they are bright enough and far enough from the quasar sight 
line). 

Given the predictions of our model, a better approach to 
studying the origin of Wo < 1 A systems at high redshifts is to 
use y-ray burst s (GRBs) as backgr ound sources rather than 
quasars [e.g., IVergani et al.l (|2009) and references therein]. 
With this technique, deep images of the fields obtained af- 
ter the GRB light has faded is more likely to reveal fainter 
galaxies, even when the GRB host galaxy is also present 
(iPrice et al.ll20MlThone et"aIll2008HChen et al.ll2009l) . 

Alternatively, a comprehensive study of absorption lines 
arising from very low redshift galaxies could be used to test 
the predictions of our model. Although HST is required to 
record spectra of UV lines at z ~ 0, nearby galaxies probed 
by quasar sightlines can be selected a priori by the same prop- 
erties studied in our model, such as their luminosity, morphol- 
ogy, star formation rate, etc. Probing a significant number of 
galaxies with a range of these properties would enable us to 
map the distribution of metal line properties directly. The in- 
stallation of COS aboard HST, with a sensitivity higher than 
the spectrographs previously available on the satellite, now 
makes a large survey possible. 

Our model gives specific predictions for the statistical re- 
lations between the absorber properties (namely, Wo) and its 
host galaxy properties (such as T and L). Nevertheless, judg- 
ing from Figure|5] the allowed parameter space in which there 
is good agreement between our model and the available Mg II 
data, is large. In this section we therefore focus on a small 
subset of all models whose parameterizations are given in ta- 
ble Q] and all provide equally good fits to the data (specifi- 
cally, all models trace the dashed black curve in Fig. |3). This 
set of models which we label A, B, and C, is consistent with 
the notion that galaxies much fainter than 0.01L* contribute 
little to the cross-section for Wo > 1A systems. This agrees 
with recent findings that galaxies much fainter than ~ 0.01L* 
are gas deficient and co ntribute little to the matter density in 
H I dZwaan et al.ll2.005b . Figure [7] shows the predicted lumi- 



nosity, star formation rate, and rest-equivalent width distri- 
butions: upper panels are relevant to situations where galax- 
ies are selected based on their associated absorber properties, 
and the lower panels to cases in which a galaxy sample is 
first obtained and follow up observations are then carried out 
to look for their associated absorption systems. As shown, 
weak absorbers that arise in disks select for fainter star form- 
ing hosts while strong absorbers tend to be associated with 
brighter star-forming objects. Nevertheless, our model does 
not preclude the possibility that some strong systems arise in 
early-type galaxies; statistically, however, the number counts 
of strong systems can be fully accounted for by star-forming 
galaxies alone. 

More quantitatively, our models suggest that the typical 
luminosities of the hosts of Wo ~ 2A systems are around 
0. 1L*. This seems t o be consistent with the recent findings 
of iChristensen et al.l (120091) who found faint L < 0. IL* galax- 
ies in the fields of such absorbers. According to the model, 
by selecting the rare, Wo > 4A absorbers, one is most likely 
to select galaxies with L> L* and with a star formation rate 
T > lM yr _1 . Nevertheless, the exact numbers depend on 
the particular model chosen, as can be seen by comparing 
models A, B, and C. For example, models B and C give very 
different predictions for the star-formation rate distribution of 
Wo > 4A hosts such that, for model C, a considerable frac- 
tion of hosts have T > 7 M yr~ 1 while, for model B, galaxies 
with T > 1 M yr~ 1 can substantially contribute to the num- 
ber of such systems. The different predictions in each case im- 
ply that reliable measurements of the L(Wo) distributions will 
constrain the model. Further, all models predict that galax- 
ies with luminosities L* <L < 6L* give rise to Wo > 4A ab- 
sorbers. The fiel ds of such ult r a-stro ng systems have recently 
been imaged by iNestor et al.l (120071) who find an excess of 
L> L* near quasar sightlines. This fact is consistent with our 
findings, although a more quantitative comparison awaits se- 
cure identifications of galaxies in those fields. 

Our model implies that surveys which first select host 
galaxies according to their luminosity and star formation 
rates, and then study absorption through their disks, would see 
a broad range of absorber properties at any given luminosity 
and SFR bin. For example, a faint dwarf with L ~ 0.01L* is 
likely to give rise to an absorber with Wo ~ 0.6A. However, 
the chances of a very bright galaxy, L > 1L* giving rise to 
the same Wo are also considerable (see Fig. |7}. This may be 
the reason for the scarcity of clear relations between galaxy 
and absorber properties (in addition to uncertainties related to 
the identification of the true host). The model is qual i tativel y 
consistent with th e rec ent findings of iBouche et al.1 (120071) . 
IWild et al.1 d2007b . and iMenard et al.1 42009J), and which are 
further discussed in §5.6. Our results highlight the need for 
large samples of quasar-galaxy pairs, which would overcome 
the large scatter in galaxy-absorber properties and help to re- 
veal the underlying galaxy-absorber association. 

The findings presented here are consistent with the notion 
that most large-scale galaxy-absorber associations, where the 
impact parameter is much larger than the H I disk size of the 
true host galaxy, arise from galaxy-galaxy clustering. This 
agrees with recent findings concerning the low covering factor 
(of order 20-50% within a few tens of kpc of the host galaxy ) 
for Mg II absorption on such scales ([Tripp & Bowed 120051: 
[ Chen & Tinker! 120081: iBarton & Cookd 120091: |Gauthier et al. 
120101 and also Bowen & Chelouche 2010). This may also be 
the reason for some edge-on galaxies apparently being asso- 




FIG. 7. — Upper panels: The predicted probability distributions for the luminosity (left) and star formation rate (right) for absorber-selected host galaxies. 
The data are shown for different Wo bins (see legend) and for the three models considered in this work (model C in red, model B in black, and model A in 
green). Clearly, a broad range of galaxy luminosities and star formation rates is selected by absorbers over a narrow Wo range, which could explain the lack of 
significant observational correlations between their properties. Generally, however, stronger systems tend to pick more luminous and more star forming galaxies. 
There is a qualitative agreement between the two models, with model B being shifted to somewhat lower galaxy luminosities and star formation rates (see text). 
Lower panels: Model predictions for the Wo-distributions as would be observed for galaxy-selected absorbers as a function of the galaxy luminosity (left) and 
star formation rate (right; note that values in the legend are in units of solar masses per year). Here too, the distributions are broad with an overall tendency for 
brighter, more star forming galaxies, to give rise to stronger systems. 



ciated with strong absorp tion along sightlin es at right-angles 
and far above their disks ( Stei del et al.l2002h . It is possible, in 
principle, to use our model and some knowledge of galaxy — 
galaxy clustering to predict the Wo distribution and covering 
factor as a function of the impact parameter from bright galax- 
ies, which are not necessarily the true hosts of th e absorber, 
and to compare to the data (IChen & Ti nker 2008, and also 
Bowen & Chelouche 2010). This raises the interesting pos- 
sibility that galaxy-galaxy clustering and group kinematics 
may be directly probed by absorber-galaxy studies leading to 
a better understanding of structure formation and to additional 
constraints on halo occupation distribution (HOD) models 
dBouche et al.ll2006t [Zheng et al.ll2007t iTinker & Chenll2008l; 
iLundgren et alj|2009t)~ 

According to the model, absorption through our own 
Galaxy disk (L ~ L* and T — lM Q yr _1 ) would yield Wo ~ 
1.2 ±0.8 A (this was calculated from the relevant Wo dis- 
tributions with the uncertainty corresponding to one stan- 
dard deviation around the average). The mean Wo or out- 
flow speed from the Galactic disk, as would be measured 
from Mg II absorption lines toward extragalactic objects, is 
~ W /2 - 0.6 ± 0.4A or ~ 60 ± 40km s" 1 (neglecting the 
possibly important contribution of galactic rotation to the ob- 
served profile along some sight lines to extragalact i c obje cts). 
This is consistent with the results of B owen et al.1 (119951 see 
their Fig. 1) who found similar velocity dispersion in several 
local Mg II systems toward quasars. We note that our analysis 
predicts that the observed velocities of the Na I lines should 
be a factor > 2 lower than those measured for Mg II (see Table 
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Model parameterization 
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(V) 


s 




fiMgll {L*) 
(kpc) 


n 


P 


color coding 


A 


0.001 


-0.7 


1.6 


40 


> i 


0.3 


green 


B 


0.01 


-0.7 


3.6 


40 


0.1 


0.3 


black 


C 


0.01 


-0.2 


2 


60 


> i 


0.42 


red 



[U; i.e., of order 30 ± 15km s _1 as would be observed along 
extragalactic sight lines. This value is q ualitatively c onsistent 
with that observed for the Ca II lines dBowenlll991l see their 
Fig. 5). A more quantitative comparison between the kine- 
matics of different ions requires better understanding of the 
structure and dynamics of galactic winds. 

Lastly, we note that our prediction for the Galactic wind 
velocity is similar to the outflow sp eed predicted at on e pres- 
sure scale-height above the disk bv lEverett et al.l (|2008) based 
on fits to the large-scale diffuse X-ray emission and using an 
energy driven outflow model. It is therefore possible that the 
outflowing cool material which is probed by Mg II lines pro- 
vides the necessary initial conditions for the onset of a large 
scale thermal wind, as advocated by these authors (see also 
§5.9). Therefore, even if intervening Mg II systems relate only 
to outflowing matter within one pressure scale-height above 
the disk, they could be related to large scale outflows that ulti- 
mately escape the galaxy and perhaps join the IGM; we return 
to this point in §5.9. 

5.2. Metallicity and Dust Content 
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FIG. 8. — The Wo — Z relation for strong absorbers at z = as predicted by 
models A (green lines), B (black curves) and C (red lines). For each model the 
mean (solid lines) and median (dashed lines) are shown. In addition, the typ- 
ical scatter in the predicted data is shown in gray points (calculated for model 
B). Although a large scatter is seen, there is a clear trend between the mean 
or median metallicity and Wo. Overlaid are data points for strong intervening 
systems for which the mean me tallicity in each bin ha s been estimated from 
the mean dust to gas ratio; see Menard & Chelouche (2009). Overall, our 
models are in good agreement with the data of Menard & Chelouche 1 2009), 
as extrapolated to low redshifts. Quantitatively, however, model predictions 
rely on the metallicity-luminosity relation whose whose applicability to ab- 
sorption systems is not well demonstrated. Furthermore, better direct esti- 
mates for the metallicity of a large sample of strong intervening systems at 
low redshifts are required make a proper comparison with model predictions. 

If strong Mg II systems arise from the ISM within, and out- 
flowing from, galactic disks, then our model can be used to 
predict the metallicity, Z, of strong absorbers, as would be ob- 
served in the local universe, using the metallicity-luminosity 
relation for z = galaxies. In particular, we take a log- 
normal distribution whose mean luminosity-dependent metal- 
licity is ~ 1.5Z (L/L*) - 5 with a scatter of order ±0.2 dex 
around the mean (Z is the solar metall i city). This is con- 
sistent with the re s ults of ITrem onti et al] d2004 see however 
iKewlev & Ellison! d2008l) for uncertainties in this relation). 
This defines the probability P(Z\L). The metallicity distri- 
bution for an absorber of a given Wo is, by construction of 
our model, P{Z\W ) = J dLP(Z\L)P{L\W ). Using these dis- 
tributions, we have carried out Monte Carlo simulations and 
constructed a mock sample of absorbers and their associated 
galaxies, which is statistically identical to that used in the 
previous sections yet now includes also the metallicity for 
each system. The results are shown in Figure [8] where a re- 
lation between the mean (or median) metallicity and Wo is 
well defined. However, the scatter is large and individual sys- 
tems having the same Wo may have their metallicities differ 
by an order of magnitude. In particular, strong systems with 
Wo > 1 A have mean metallicities of (Z) > Z Q (e.g., Wo ~ 5A 
systems have (Z) > 2Z ) while weaker systems, which origi- 
nate in disks, have sub-solar metallicities. Evidently, the exact 
numbers depend on the particular model chosen. We empha- 
size that these results are model predictions for the metallicity 
of strong Mg II systems, as would be observed at z = 0. Pre- 
dictions for higher redshift systems will require knowledge of 
the metallicity of high redshift disks. A quantitative extrapo- 
lation of our results to high redshifts is beyond the scope of 
the present work. 

Comparing to data in the literature, we assume that the dust- 
to-gas ratio of strong intervening systems can be used as a 
proxy for their metallicity. We therefore take the statistically 
measured dust-to-H I column r atio, as converted to metallic - 
ity and extrapolated to z — by Menard & Chelouche (2009), 



and plot the mean ratio and its uncertainty in Figure [8] As 
can be seen in the figure, both the mean and median metal- 
licities derived using model C predict values that are larger 
than observed. Models A and B, however, agree with current 
measurements surprisingly well. In particular, the predicted 
trend whereby weaker systems have a lower mean metallicity 
is in nice agreement with the observations, and the predicted 
values are consistent (up to ~ 20% for the mean) with the 
data within the quoted uncertainties. The obtained Wo — (Z) 
relation may also explain the relatively rapid increase in the 
equivalent widths of Fe II A 260 absorption line with Wo Tat a 
given H I column; see Fig. 1 in Menard & Chelouche (2009)] 
since stronger systems are more metal-rich. 

The agreement between model predictions and the extrap- 
olated data is surprising given the crude metallicity estima- 
tor s used. For ex a mple, the metallicity-luminosity relation 
of ITremonti et al.l ((2004) was determined from metallicity 
measurements of H II regions on compact scales in galax- 
ies. However, it is not clear whether the metallicities of 
gas in galactic H II regions need be the same as interstel- 
lar gas between H II regions . In t his respect, it is interest- 
ing to note that iBowen et al.l (120051) found good agreement 
between the metallicity of a Damped Lya absorption (DLA) 
system [absorpt ion by neutral hydr ogen with column densities 
logA^Hi > 20.3 dWolfe et al.lll986h l measured from its metal 
absorption spectrum and that deduced from emission. There 
are, however, contrary examples where H II region metal- 
licities do not agree with interstellar absorption metallicities 
(ICannon et all 120051: Icliurchill et all 120051) . alt hough untan- 
gling the effects of metal licity radial profiles dGiveon et al.l 
l2002t iBresolin et al.l 2004) — and indeed, deciding whether 
the correct galaxy has been identified as the origin of the DLA 
system — complicates the task of matching emission and 
absorption metallicities. Our models do not include metal- 
licity gradients, an omission that will lower the predicted 
metallicities for sight lines that intercept the outer regions of 
disks. While metallicity gradients can have a crucial effect 
on the large-scale extinction through extended galaxy disks 
(lZaritskvlll994t iHolwerda et al.ll2005l IChelouche et alJl2057t 
iMenard et al.ll2010l) and may be easily incorporated into our 
calculation scheme, we believe that more elaborate models 
are currently unwarranted as the involved uncert ainties on the 
metal l icity of disks on la rge scales are large dGiveon et alj 
120021; iBoissier et al.l 120041) . Better estimations of the metal- 
licity or dust-to-gas ratio of Wo > 0.5A systems promise to 
shed light on the metallicity gradients in the outer regions of 
disks. 

DLAs are known to c onsistently show me t allicities less 
than the solar value r e.g.. | Petitiean et all d2000l) ; iNestor et al. 
(2003); iKhare et alj d2004l) ; IChenTtTiT(l2005l) ; iPeroux et al. 
(2007)]. While at first this seems to contradict our model pre- 
dictions, we note that the probability for obtaining a DLA for 
a strong intervening system with a given Wo is similar, within 
a narrow Wo range, for all systems with Wo > 0.6A dRao et al.l 
2006) with no DLAs currently selected by Wo < 0.6A sys- 
tems. Therefore, most DLA surveys, which select by H I col- 
umn density may, in fact, select 0.6 < Wo < 1 .5A systems be- 
ing the most abundant (see Fig. [3]). The metallicity associated 
with such systems, in our model, is indeed sub-solar. Further- 
more, DLA systems are observed at high redshifts while our 
formalism was developed (and applied for extrapolated Mg II 
data) at z — 0. If the evolution in metallicity of DLA systems 
is similar to that of Mg II systems then we expect their typical 
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FIG. 9. — The number density of Mg II systems per unit redshift, dN(z = 
0)/dz with a H I column density above a threshold, A^/j™, as expected for 
strong systems with Wo > 0.5A (black squares). This was calculated using 
Q l dW (d 2 N/dz3Wo)P Wo (N m > Aftf where P Wo (N Hl > A$f ) is the prob- 
ability for an Mg II absorber with rest equivalent width Wo to have a H 1 col- 
u mn density greater than A/jjj 1 . This probability was estimated from the data 
of Raoetal. 12006, see their Fig. 2) and is shown in the inset. Also shown 
in dashed line are the dN(z = Q)/dz as measured for nearby galaxy disks 
bv lZwaa n et al. 12005). The uncertaint y was taken to reflect p otentially im- 
portant systematic effects discussed in Rvan-Weber et al. (2003) and shown 
as a gray area. Evidently, galaxy disks at z = can fully account for the 
(expected) number of z, = Mg II systems up to the highest column density 
bin. However, the number counts of strong systems with the highest columns, 
N m > 10 213 cm 2 , is expected to be a factor ~ 2 — 3 higher than estimated 
based on 2 1 cm emission line observation of galaxy disks (see text). That 
said, these estimations are sensitive to beam size effects as well as extinction 
bias and evolution. A crude correction for the latter effects may go either way, 
as shown by the arrows (potentially important correction factors for the lower 
column points are not shown). Error bars reflect the uncertainty on the pre- 
dicte d absorber signa l due to Poisson errors in determining P(Nm > A'Jjj"; Wo) 
from Rao et al. 1 2006) data set and does not include additional uncertainties 
related to the form of d 2 N(z = 0) /dzdWo- 

metallicity to be lower still . Specifically, using the evolution 
in the dust-to-gas ratio from Menar d et al.1 (120081 see their Fig. 
3) as a proxy to the evolution in the metallicity, one finds that 
typical DLA metallicity should be < 0.2Z e at z > 1.5. At 
still higher re dshifts z ~ 3, and usin g the metallicity evolu- 
tion curves of iKulkarni & Fali (120021) . the typical metallicity 
is exp ected to be < 0.02Zp , and in agreement with observa- 
tions (IKulkarni et al.ll2005h . We note, however, that the scat- 
ter around the mean value is expected to be large (see Fig. [SJ, 
which cou ld account for the large r ange of metallicities re - 
por ted bvlPeroux et al.l (120071 120081) . iFrank & Perouxl (12010b . 
and Kapla n et al.1 (120101) . In addition, extinction effects could 
bias magnitude-limited samples agai nst finding hig h metallic- 
ity, high dust content DLA systems (Bouche 2008, and refer- 
ences therein). 

5.3. The H I distribution of strong Mg II Systems 

If disk galaxies give rise to Mg II absorption, then the ob- 
served Nm distribution of strong systems should be accounted 
for by disks. At present, our knowledge of the H I column 
densities associ ated with stron g Mg II systems is limited to 
high redshifts dRao et al.ll2006l) . while our knowledge of the 
H I distribution in disks is limited to nearby objects. Here 
we make use of the probability for an absorber of a given 
Wq to show a HI column density greater than some value, 
Pw (N m > N^). This probability was determined from the 
observations for z ~ 1 systems by iRao et all (120061 see their 



Fig. 2) and is shown in the inset of Figure @H 

Here we assume that the above determined Pw (Nm > A^J 1 ) 
is applicable to low-z systems and calculate the 
number density of Mg II systems per unit redshift 
that give rise to N m > A$f, dN(N m > A$f)/dz = 
Q A dW (d 2 N/dzdWo)P Wo (N m > N 1 ™). We take 

d 2 N /dz dWp from our be st-fit model to the extrapolated 
data of dNestor et al. 2005) and show the results in Figure 
|9] The uncertainties on the predicted dN(Nm > N^)/dz 
signal shown in Figure [9] include only those due to Poisson 
noise in determining P\y (Nm > ^h" 1 ) fr° m the IRao et all 
(2006) data and do not include additional uncertainties of a 
factor > 3 on the amplitude of the d 2 N /dzdWp signal (see 
§3). Comparison with the H I data of IZwaan et al.1 ( 2005, 
see their Fig. 8), which were obtained from 21 cm emission 
observations of local galaxy disks shows overall agreement: 
the shape of the dN(Nm > N m r ')/dz curve is reproduced 
given the uncertainties. Furthermore, the amplitude is also 
consistent within the uncertainties lending further credibility 
to the model and our estimates for the cross-section of 
galaxy disks. The uncertainties on the 21 cm e mission data 
show both the statistical errors reported by IZwaan et al.l 
(120051) as well as possible systematic effects discussed 
in iRvan- Weber etal] d200l . For JV|f < 10 19 8 cm- 2 , the 
WHISP sample suffer from sensitivity issues and is therefore 
incomplete. However, an extrapolation of the powerlaw trend 
at low column densities as implied by 21 cm observations is 
seen to be consistent with the much lower columns probed 
by strong Mg II systems. This may be related to the fact that 
strong systems seem to trace the Galactic dust-t o-gas relation 
( Bohl in et a!J1978t iMenard & Chelouchdl2009h . 

As is evident from Figure [9] the number counts of strong 
systems with the highest Nm seem to be higher by a factor 
of ~ 3 than implied by the 2 1 cm emis sion data. This fact 
was already noted by IZwaan et al.l (120051) and it remains to 
be seen whether it is persists when better statistics are gath- 
ered with COS. The reason for this mismatch is not known 
and several possible factors that could cause it: for exam- 
ple, if the local Schmidt-Kennicutt law holds also at high 
redshift then the fact that the global star formation rate in - 
creases with redshift dLillv etal.lll996t iMadau et all 119981) . 
may mean that Nm increases with redshift, possibly in a com- 
plicated manner. Therefore, it may be that our assumption 
that Pw {Nm > A^j™), as derived for high-z systems, is di- 
rectly applicable to the local universe is unwarranted. Specif- 
ically, given the fact that the global star formation rate has de- 
creased by about an order of magnitude from z ~ 1 to present 
times, and using the Schmidt-Kennicutt scaling A^hi « Y 2 / 3 , 
A^hi may be lower by a factor ~ 4 compared to that mea- 
sured by IRao et al.l (120061 see Fig. |9]l for high-z systems. A 
low-z analog to the IRao et all (120061) compilation or higher- 
z 21 cm observations of galaxy disks will be most useful in 
settling this issue. A further systematic effect may be re- 
lated to the typical beam size used by each approach: in ab- 
sorption line studies the quasar emission region sets the an- 
gular resolution, which is of order ~ 10~ 8 rad in diameter. 
In contrast, radio surveys have typical beam sizes which are 
~ 10 3 times larger. As demonstrated by IZwaan et al.1 (120051 

8 We note that while Rao et al. 1 2006) use various selection criteria to con- 
struct their sample, which particular set of criteria used is not important and 
results in a similar Pw (Nm > ^nC) given the uncertainties. 
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and references therein), lower resolution results in an under- 
estimation (smearing) of high column density regions. The 
degree to which this explains the apparent mismatch in the 
number of systems with the highest H I columns is unclear 
and depends on the structure of the ISM. Lastly, it may also 
be that th e high est Nya systems are over-represented in the 
iRao et all (12006b sample due to the combination of selection 
criteria used; an effect which cannot be easily identified given 
the small number statistics for jVhi > 10 21 cm~ 2 systems in 
their work. We note that a recent census of DLA systems by 
Noterdae me et al.1 (l2009|) at z > 2 shows c onsistency, up to a 
scaling factor, with the lZwaan et al.l d2005) results. 

It is worth noting that there is an additional, potentially 
important, syst ematic effect t hat re lates to dust extinction 
and reddening: Me nard et al.1 (120081) have shown that high- 
column systems may be missed by magnitude-limited sam- 
ples, such as the Sloan Digital Sky Survey (SDSS). Us- 
ing the mean dust-t o -gas ratio for strong systems from 
iMenard & Chelouchel (120091) and taking an SMC reddening 
curve, we estimate E(B — V) ~ 0.2 mag for Nm ~ 10 213 cm~ 2 
systems, implying a mere ~ 50% comp leteness of the SDSS 
for such systems (Me nard et al.l 120081 see their Fig. 2). 
Higher column density regions, assuming the same mixture 
of dust and gas, may be completely missed by magnitude- 
limited samples, leading to further uncertainties concerning 
the true column density distribution of strong systems. It is 
important to try to map the outer regi ons of galaxies in 21 cm 
down to the lowest possible columns dBraun & Thilkerl l2004l) 
and test whether the cross-section and column density distri- 
bution of disks and their outflows at low-z can indeed explain 
the properties of strong systems. 

Our finding that most strong Mg II systems are account- 
able for by disks implies that one does not lose most strong 
absorption systems due to an extinction bias, as concluded by 
Ellis on et all (12004^1 In our model, the DLA and strong Mg II 
systems phenomena are consistent with being two facets of 
the same phenomenon: ISM in (and outflowing from) galaxy 
disks. This is in line with the results of other studies that com- 
pared the properties of H I and Mg II selected DLAs finding 
no statistically significant difference between the two popu- 
lations (Pe roux et al.ll2004h . It is also in agreement with the 
cosmic mass density of Mg II selected DL As and H I selected 
DLAs being similar [see e.g., Fig. 17 of iNoterdaeme et aTl 
(2009)]. Our findings are also consistent with the notion that 
the high-end tail of the observed H I column density distri- 
bution, for a given Wo, arises in le ss luminous, hence metal- 
and dust-poor disks (Bouche 2008) while sight lines through 
metal-rich star forming galaxies having, on average, large H I 
columns and dust-rich gas, may be missed due to extinction. 
If true then, on average, galaxy disks allow through a consid- 
erable fraction of the the light from background sources. This 
conclusion can be checked, observationally, by deriving a bet- 
ter estimate for the local value of d 2 N /dzdWo (especially at 
the high-Wo end of the distribution which would be most af- 
fected by the extinction bias), and by carefully searching for 
background sources behind crowded stellar fields so that their 
density may be compared to that of objects in the field. A 
better understanding of the opacity of star-forming disks to 
background radiation would also reveal the degree to which 
our parameter estimates are affected by extinction bias at the 

9 However, one may still lose a large fraction or all of the less numerous 
ultra-strong systems or those fewer systems with large H I columns; see §5.2 
and also Menard et aT H2008l) . 



high-Wb and high-T tails of the distributions, and would help 
to determine whether additional physics (such as a luminosity 
dependent /l,^w) is required. 

5.4. Starburst Winds and Mg II Systems ' Kinematics 

Our model implies that the number counts of strong in- 
tervening systems can be fully accounted for by the num- 
be r of star forming di sks [/z, ~ 1 given the number counts 
of iMarzke et al.l (1 1998b 1 if the cross-section for strong Mg II 
systems is of order the area of the H I disk ~ 1). This 
allows us to calculate the mean Mg II absorption line pro- 
file that would be observed against the galactic continuum 
for an ensemble of galaxies (e.g., with a given star formation 
rate) by stacking their individual (predicted) spectra. Specifi- 
cally, using the luminosity function of galaxies and their star- 
formation probability distributions, we created a mock cat- 
alog of galaxies and their associated absorbers according to 
the modeling scheme of §3. An absorption spectrum for each 
galaxy was then constructed under the assumption that the ab- 
sorption trough is saturated and t hat the outflow fully co v- 
ers the UV emitting disk [but see lMartin & Bouchel ([2009)1. 
This means that we approximate the absorption line profile, 
as would be seen along a single sight line through a disk, as a 
rectangle centered at the systemic velocity of the galaxy. The 
mean line profile for a sample of galaxies, with some given 
properties, was then calculated by taking the mean flux at 
every photon energy of all spectra. Th is procedure is qual- 
itatively similar to that carried out by iWeiner et al.l d2009l) 
for DE EP2 data. To allow comparison with the Weiner et al. 
(2009) results, we calculated v'75%, the velocity at which 75% 
of the flux is transmitted through the blue wing of the (mean) 
line profile, as a function of the SFR of the galaxies . 

The resulting line profiles are shown in Figure [10] The 
shape of the predicted line profile, as would be observed along 
sight lines to background objects that pierce through interven- 
ing disks, is characterized by a deep trough at the systemic ve- 
locity of the galaxy with symmetrically extended blue and red 
wings. The shape of the mean line profile is not determined 
by optical depth effects — individual systems are saturated by 
construction — but rather by effective partial covering effects 
when averaging over the full sample of galaxies. This comes 
about from the fact that some systems in the ensemble have 
high velocity absorption components while others do not. In 
the latter case the unattenuated continuum contributes to the 
mean line profile. This predicts saturated mean absorption 
line troughs that do not reach a flux level of zero, consis tent 
with the 1 : 1 doublet ratio reported by Weine r et al.l d2009l) for 
their mean line profiles (note that our modeled absorption line 
profiles are not convolved with the instrumental line spread 
function hence the black troughs at the systemic velocity). 

The exact velocity up to which the line profile extends is 
not very sensitive to the chosen parameterization (solid line 
models in Fig. flOt . Nevertheless, for one of our runs, we as- 
sumed that the UV emitting region is of radius po = rjRMgiu 
i.e., smaller than the size of the H I disk. In this case only 
the outflow kinematics within po will contribute to the for- 
mation of the mean line profile (magenta line in Fig. ITOb . 
Models in which the outflow kinematics varies across the disk 
(e.g., model B) result in the absorption trough extending to 
much higher velocities. More quantitatively, for model B, we 
find that, for 10 < T < 13M e yr" 1 (30 < T < 6OM yr" 1 ), 
V75% ~ 300km s (v75% — 450 km s ), which is consistent 
with the results of I Weiner et al.l ([2009) for Mg II winds (see 
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FIG. 10. — Mean Mg II A2796 absorption line profiles from galaxies with a 
specified range of star formation rates (see legend). Color coding of the dif- 
ferent models follows that given in table[T] In addition to models A, B, and 
C, we consider a model for which the UV emitting disk size corresponds to 
the core radius while the outer parts of the disk are UV faint and do not con- 
tribute to shaping the mean profile (magenta lines). The trough is symmetric 
but only the blue wing of the line is shown. The line shape is determined by an 
effective partial covering effect resulting from averaging over many systems, 
some of which having broader troughs than others (see text). Model B results 
in broader troughs on accounts of its higher S, v value compared to model C. 
The results for the velocity at 75% transmission (see horizontal dotted line) 
for model B, with its core radius being the size of the UV emitting disk, are 
in agreement with the measurements of Weiner et al. 1 2009) for galaxies with 
10 < X < 6OM yr~ ! whose mean v'75% and vgo% are shown as black points 
and full velocity range shown with error bars. 

Fig. fTob . Results for models A, B, and C for which the UV 
emitting surface is the entire disk or the wind velocity is con- 
stant across the disk, seem to under-predict the observed char- 
acteristic outflow velocities. Our results imply that galaxies 
with a relatively l ow SFR, will sho w only low to moderate 
velocity outflows (Rubin et al. 2009)). 

When observed at high spectral resolution, strong interven- 
ing systems often show several discrete absorption compo- 
nents, implying that the absorbing medium is clumpy. Draw- 
ing from the apparent correspondence between intervening 
Mg II systems and galactic winds, this may indicate that galac- 
tic winds are largely clumpy and that their structure is char- 
acterized by clouds and filaments. Some support for this no- 
tion COTn^sJrcm NcH s |u^^ 



ies (Rupke et al 


2002; Martin 2005:lMartin & Bouchel!2009; 


iSato etalJ 12009 


). Indeed, a highly complex density struc- 



ture may be expected in such dense and dynamically evolv- 
ing environments due to the development of various insta- 
bilities ( e.g., Kelvin-Helmholtz) as the winds accelerate an d 
expand dStrickland & Stevensl [19991: IStrickland etail 120041) . 
Conversely, the physical picture of galactic outflows naturally 
explains the rich line profile phenomenology of intervening 
systems: in this picture, the rich multi-component absorp- 
tion that characterizes intervening systems is mostly due to 
(hydro-) dynamical effects related to outflow physics and is 
not related to the clustering of halos. One way to test and 
refine our model would be compare the absorption line pro- 
files of galactic winds, as observed toward the centers of star- 
forming disks, with those seen toward background objects 
whose sight lines pierce through the disks of the same galax- 
ies. Even more revealing would be to find background objects 
behind star-forming galaxies and use their spectra to constrain 
the small scale physics of galactic winds, which is not acces- 
sible for studies using the galaxy itself as a background o n 
accounts of the large emitting area (Martin & Bouche 2009). 
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FIG. 11. — Evolution in the statistical properties of galaxies with cosmic 
time. Shown in solid (dashed) line is L*(z)/L*(z = 0) ((j>* (z) / <j>* (z = 0)) 
ba sed on a second orde r polynomial fit to t he L (z) (<j>*(z)) m easurements 
of IGabasch eTaTI J2004l) (for z > 0.45) and IJames efaH 12503) (for z ~ 0) 
in the 5-band. If the luminosity of galaxies at any epoch is a measure of 
the star formation rate then the evolution in the star-formation rate density, 
Psfr (z) I Psfr (z = 0) is consistent with the observations. D ata points for 
Psfr (z) / Psfr (z = 0) at z < 3 were taken fro mlLv et all f2 007 ) and were nor- 
malized according to the measurements of James et al. (2008) at t he local 
univer se (filled circles). High redshift data w ere taken from IHeayensetaT] 
1 2004) based on data from Steidel et al. 1 1996, empty circles) and Ouchi et al. 
1 2004, empty squares). Unsurprisingly, our parameterization for the luminos- 
ity function traces the global star formation rate density. 

5.5. Redshift Evolution 

Several authors have found indications for evolution in the 
numb er density of intervening systems as a function of red- 
shift dSargent et al.ll 19881: iPetitiean & Bergeronlfl990t ISteidell 
1990). Mos t relev ant to our study are the recent results of 
Nesto r et al.1 (120051) who find tentative evidence for redshift 
evolution in the number density of absorbers, d 2 N/dzdWo, 
as a function of Wo. In particular, evolution effects seem to 
be more pronounced for strong absorbers and the number of 
such systems decreases from high r edshifts to present times 
[see, however, Lun dgren et al.l (120091) who find that the num- 
ber density of Wo > 2A systems actually increases from high 
redshift to z — 0.4]. Nevertheless, the uncertainties in the ob- 
served evolution of d 2 N / dzdWo are large and better statistics, 
as a function of Wo, are required especially at low redshifts. 
In our model, Wo is a direct consequence of the SFR of galax- 
ies, which is known to be evolving. It is therefore natural 
to ask what our model predicts for the redshift evolution of 
d 2 N/dzdWo given the evolution in the properties of galaxies 
over cosmic time. 

In its present form, our model does not include evolutionary 
effects and additional assumptions are required to incorporate 
them into the formalism if any comparison is to be made with 
data relevant to absorption line systems at high-z. As noted 
in §3, the properties of low luminosity star-forming galaxies 
are poorly determined at high redshifts. Nevertheless, using a 
specific set of assumptions we can check whether our model 
is consistent w i th obs ervations. We rely on the recent work of 
Gabas cTTet alj ([2004 ) who find that the Schechter luminosity 
function provides an adequate description of galaxy statistics 
at high-z, and give values for (j>* and L* for several redshift 
intervals. Due to the slope of the faint end of the high-z lu- 
minosity function being difficult to determine due to incom- 
pleteness, we hold a fixed at the local value. (A similar ap- 
proach was taken by Gabasch et al. only with a = —1.25.) 
We fit (j)*(z) and L*(z) values reported b y Gabasch et al. for 
0.45 < z < 5 together with the results of Uames et alj (120081) 
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FIG. 12. — Evolution in the statistical properties of strong intervening s ystems. Left: The dN/dz evolution with redshift for several values of Wo > W(J lm as 
denoted next to each curve. Also shown are data from Nestor et al. 12005) for strong systems in the redshift range 0.4 < z < 2.3. Two versions of model B are 
considered (solid and dashed lines) that are both consistent with the shape and amplitude of the signal for Wo > lA systems. In particular, all models predict 
a sharp rise in Wo > lA systems' density toward z ~ 2 followed by a gradual fall which follows from the decline in the global star-formation rate. Stronger 
systems number density peaks at higher redshifts. More specifically, our model B parameterization with £ w = 3.6, RMgll = 40kpc (dashed line) over-predicts 
the signal for the strongest systems by 15-50%. A better fit to current data over the full Wo and redshift range is obtained for = 3.0, Ruigll — 43 kpc. The 
model predicts very different evolution track for strong systems than non-evolutionary models (dashed blue line for the case of Wo > 1 .5A systems. High redshift 
data are scarce and our model is consistent with current data [a point from Jiang et al. (2007) for z ~ 4.3 systems is shown]. Right: Our modified model A 
(dashed line) and model B (solid line) predictions for d'N(Wo,z) /dzdWo for several redshift epochs (see color coding) shows evolution which is prominent for 
the strongest systems. In particular, the distributions functions flatten at high redshift showing the abundance of very strong systems at high redshifts. At low 
redshifts, the contribution of weaker systems (with Wo ^ lA) is relatively large while that of super-strong systems (Wo ~ 4A) drops considerably, by about an 
order of magnitude, from z = 0.5 to present times. Also shown is a departure from a pure powerlaw at the high-Wo end of the distribution for our revised model 
B (as compared to our revised model A; dash dotted line), resulting from the contribution of sight lines passing close to the centers of relatively low luminosity 
galaxies. Testing which of the models is relevant requires better knowledge of the kinematics of galactic outflows across disks and better statistics for strong 
systems. 



for z = by a second order polynomial so that: 

log</>*(z) = -2.65[log(l+z)] 2 +1.031og(l+z)-2.52 
logL*(z) = -0.30[log(l+z)] 2 + 1.421og(l+z) + 10.14 

(11) 

These relations, normalized to local values, are shown in Fig- 
ure HU 

Incorporating redshift evolution into our model is done 
in a straightforward way: we assume that all local rela- 
tions and distributions hold at high redshifts and that only 
the luminosity function changes with cosmic time via equa- 
tion ITTb . Some support for our simplified approach may 
be gained from the findings of ( Ravindranath et al.l 120041) 
who found that the luminosity-dependent cross-section (or ra- 
dius) of galaxies does not evolve with redshift in the range 
0. 25 < z < 1, in agreeme nt with the theoretical expectations 
of lFall & Efstathioul(ll980t) . 

If the characteristic SFR of an L* gala xy is identified, o n 
average, with its luminosity, at any epoch dJames et al.ll2008l) . 
then the predicted global star formation rate density of the 
universe, Psfr> based on our parameterization for the lumi- 
nosity function evol ution with red shift, is consistent with cur- 
rent measurements dLv et al.l2007l and references therein; see 
Fig. ITTb . We caution, however, that some of the SFR measure- 
ments have recently been questioned and variou s systematic 
effect s may exist in the estimation of Psfr(^) dNordon et alj 
2010). The scenario defined above directly links the evolu- 
tion in the statistical properties of strong intervening systems 
to that of star-forming galaxies and is therefore quali t atively 
different from the one advocated by Tinker & Chenl (120091) 
in which absorbers' evolution is linked to the growth of dark 
matter halos over cosmic time. 

We integrated equation|5]to calculate dN(Wo > WQ lm ;z)/dz 



for our model B for several values of W llm — see Figure [T2] 
Our predictions for dN/dz are very different from the ex- 
pectations of non-evolutionary curves [« (l+z) 2 /£(z)]. In 
particular, the decline in the number density of absorbers to- 
ward low redshifts is more pronounced than the one predicted 
by non-evolutionary curves by a factor of a few. This fol- 
lows from the rapid decline in the star-formation rate from 
its peak value at z > 2. At still higher redshifts, the mean 
star-formation rate density declines as does the predicted 
dN(z) /dz, which traces it, and is therefore very different than 
the predictions of the non-evolutionary curves (see Fig. [T2l . 
It is worth noting that a very different model, which ties the 
evolution of strong intervening sys tems to that of intermedi- 
ate mass halos, was sugg ested by iTinker & Chenl (120091 see 
their Fig. 6) and predicts a much steeper decline in dN/dz 
for z > 2 than the one found here. Our calculations show 
that the density of absorbers above some Wo limit peaks at 
intermediate redshifts, depending on Wq (Fig. [T2l . Specif- 
ically, the density of Wo > lA systems peaks at z > 2 while 
that of stronger systems with Wq > 3. 5 A peaks at earlier times 
corresponding to z > 3.5. At present, high-z data are scarce 
and our model is consistent wit h recent measurem ents of the 
dN(W > 1.5A:z~4.3)/dz bv l3lang^n(l2007h . 

Comparing our model predictions to observations we see 
that the curvature in the data as a function of redshift is well 
traced by the models for a broad range of W llm . Furthermore, 
while the amplitude of the signal for Wo > 1 .5 A systems is 
matched by the model, the model seems to over-predict the 
signal for W > 2. 5 A by ~ 15% an d that for W > 3. 5 A by 
> 50%. As noted by iNestor et all d2005l) . the statistics of 
Wo > 3. 5 A systems is poor and quantifying the systematics 
that affect the dN/dz measurement is challenging. It is there- 
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fore possible that the true number density of such systems is 
actually higher than reported in their work. Whether this can 
account for the entire difference between the dN / dz predicted 
by model B and the observations is unknown. Alternatively, 
our model, whose normalization relies on extrapolated data 
from high-z measurements, may be inaccurate. As discussed 
in §3, we have ext rapolated the non-evolutionary curves of 
iNestor et al.l (|2005) to z = to determine the relative normal- 
ization of the d 2 N(W ,z = 0)/dzdW Q curve while our model 
predicts a more rapid decline in the number of s t ronger sys- 
tems at low-z. Attempting to fit the INestor et al.l (120051) data 
set for dN(Wo > W llm ;z)/<9z, we arrive to a somewhat dif- 
ferent parameterization than our model B, which we term as 
model Bl, with % w = 3, R M gn(L*(z = 0)) = 43 kpc and with 
all other parameters held fixed at their model B values. This 
model shown in solid line traces well the shape and amplitude 
of the observed signal for a broad range of W lim . (We note 
that the model is also in good agreement with the other data 
sets considered in the previous sections.) Using model Bl to 
match the dN(W > W^ m ;z)/dz data of INestor et ail (120051) 
over a broad redshift range provides us with a refined estimate 
for d 2 N(Wo,z — 0)/ dzdWo which is characterized by a lower 
density of the strongest systems that suffer from the mostpro- 
nounced evolutionary effects at low redshifts; see figure [5] 

As already noted in §3, predictions for 
d 2 N(Wo,z = 0)/dzdWo are expected to show deviations 
from a pure exponential for models of type B in which rj < 1 . 
Studying the shape of d 2 N(Wo 7 z)/dzdWo as a function 
of redshift shows that this is indeed the case and that the 
contribution of relatively faint galaxies to the high Wo-end 
of the distribution is most pronounced at z = 0, occurring 
for Wo > 4A systems where the contribution of the relatively 
rare star-forming systems is small. At higher redshifts the 
deviation from a powerlaw is less pronounced yet may be 
possible to detect even at z ~ 1 provided the statistics of 
Wo > 5 A systems may be determined to better than ~ 10%. 
For comparison we also plot a version of our model A with 
E, w = 1.37 and all other parameters fixed at their model 
A values. This model provides a good fit to the d ata sets 
described above (aside from the ( Wein eret al.ll2009l) results; 
§5.4) but does not show deviation from an exponential decay 
in d 2 N(Wo,z)/dzdWo even at large Wq. Therefore, a test for 
the connection between outflows from star-forming disks and 
strong intervening systems would be to check for deviations 
from an exponential at high Wo values at high redshifts 
and compare those to the mean Mg II absorption troughs 
observed toward star-forming galaxies. Nevertheless, any 
proper comparison between the data and the model requires 
a good understanding of potential systematics, such as dust 
extinction, that are expected to be more important at small 
impact parameters from galaxy centers, and which would 
select against detecting high Wo systems. 

If this model is correct then by measuring the density of 
absorbers over a broad redshift range one may be able to 
measure Psfr m a wa Y which is complementary to that ac- 
complished by studying the emission characteristics of star- 
forming regions. In fact, better measurement of dN(Wo > 
W llm ;z ~ 4)/<9z may provide interesting constraints on the 
SFR history of the universe and help constrain the physical 
model proposed here. It is, however, important to stress that, 
in order to be able to provide robust star formation density 
measurements, the physics of galaxies and their associated 



outflows needs to be better understood. Alternatively, should 
Psfr be estimated by other means, it may be possible to inves- 
tigate the physics of bright as well as faint galaxies, in terms 
of their outflows and disk sizes, up to the highest possible 
redshifts (so long as a bright background object is present) 
thereby evading the limitations plaguing emission-based stud- 
ies. 

5.6. Star Formation Rates for Strong Systems 

Several studies have aimed to measure the star forma- 
tion rate associated with intervening absorption line systems 
dWolfe et al.ll200l|Kulkarni et al.ll2006tlWorfe & Chenl | 2006 | : 
[ Bouche et all 120071: iGharanfoli et all 120071: iWild et al J 120071: 
Noterdaem e et al.H2.Q10t INestor et al.H2.Q10t) . In some studies, 
only upper limits for the star formation rates of host galaxies 
were obtained, while in other, lower bounds were measured 
due to contamination by the quasar light, finite fiber or slit 
sizes, and dust extinction (assuming that the true host galaxy 
had been detected). In a few cases, integral field unit (IFU) 
observations were obtained leading to more reliable estimates 
for the SFR associated with strong systems. 

We have used our model Bl at various redshifts, as de- 
fined above, to calculate the SFR distributions as a function of 
Wo. In particular, the results for a run of our model with 10 4 
absorber-galaxy pairs are shown in Figure [T3] where although 
most absorbing systems are associated with low star-forming 
galaxies, the scatter is large and a high star-formation tail is 
observed. This behavior is reflected in the large difference 
between the median and the mean of the SFR (as calculated 
from a run of our model Bl with 10 6 absorber-galaxy sys- 
tems) especially for Wo > 1 A absorbers. A clear trend is seen 
whereby the mean or median SFR increases with Wq. The 
mean and median values are also plotted for higher redshift 
models. These SFRs show a slight increase with z. Below, 
we discuss some of the datasets used to derive SFRs at differ- 
ent redshifts, and examine how consistent they are with our 
models. 

Measurements of the star formation rate for i ndividual sys- 
tems w ith known Wq have been reported by Bouche et al. 
d2007b . using IFU observations, and are shown in Figure [T3l 
These data correspond to z — 1 systems with which our model 
predictions are consistent . Additional meas u remen ts have re- 
cently been reported by iNoterdaeme et"ail (2010) for lower 
redshift (z ~ 0.5) systems and our model predicts compara- 
ble mean and scatter around the mean (compare the scatter in 
the simulated systems to the measurements in Figure [PI ). It 
should be noted that while the INoterdaeme et"ai] (|2010) study 
is biased to galaxies with brighter emission lines, the reported 
SFR should still be treated as a lowe r limit on accounts of fiber 
losses. Just how representative the Noterda eme et alj (120101) 
measurements are of the true SFR in absorber hosts is yet to be 
determined. Very high star formation rates around two sight 
lines with ul tra-strong sy s tems a t z ~ 0.7 with Wo > 3A were 
reported by INestor et al.l d2010l) . consistent with our model 
predictions. We note, however, that in these cases an interact- 
ing galaxy pair is observed and it is not clear which galaxy 
hosts the absorber and whether a particular galaxy-absorber 
a ssociation is at all mea ningful in these cases. 

Kulkarniet al. (2006) placed upper limits on the star for- 
mation rate in z ~ 2.4 DLA systems of > 1 M© yr . As dis- 
cussed in §5.3, a large fraction of the DLA population may 
be associated with strong Mg II systems. Given the results 
of Ra o et al] ([2006) that show a rather uniform column den- 
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FIG. 13. — Model predictions for the mean (solid line) and median (dashed 
line) star formation rate for z = (black curve), z = 1 (red curve), and z = 2 
(dark red curve) models. For the z = model, a plot of the simulated systems 
is shown exemplifying the large scatter and the high SFR tail of the distribu- 
tion resulting in the mean values being considerably higher than the median 
values, especial ly for low Wp . Also shown are measurements for individual 
absorbers from Bouche et al. 12007) whose average redshift is unity (filled 
circles and filled upper limit triangles). Empty triangles represent upper limit 
data for z > 2 DLA systems from Kulkarni et al. 12006), for which Wo is not 
known and a value of Wo = 0.6A is assigned (see text). The uncertainty on Wo 
in this case was calculated from the d 2 N(Wo,z = 2)/dzdWo distribution and 
marks the Wo range in which 68% of the syst ems are expected to lie. Recent 
measurements from Noterdaeme et al. (2010) are shown in gray squares. We 
also show data from Wild et al. |2007) which correspond to a statistical mea- 
surement at low (0.4 < z < 0.8; small star symbol) and high (0.8 < z < 1.3; 
large star symbol) redshifts, and for which we also assign a value of Wo = 0.6 
and calculate the uncertainty as before. New dust-corrected measurements 
for the star-formati on in the environment of ultra-strong z ~ 0.7 systems by 
INestoretaU l2010) are marked as stars. Evidently, our model predictions are 
consistent with the data for the relevant redshift range. 

sity distribution for systems of dif f erent W p (for Wo > 0.6A), 
the objects in the Kulk arni et alJ (l2006h sample are likely 
to be associated with the more numerous Wo > 0.6A sys- 
tems [for W) < 0.6A systems the probability for detecting 
a DLA is consi derably reduced given the data reported by 
Rao et al. (2006)]. Figure[13]compares the upper limits of the 
Kulkarni et alJ (2006) sample, to which we have assigned a 
value of W) = 0.6A (and note that 68% of selected systems are 
expected to lie in the range 0.6<W < 1.7A given the shape 
of the d 2 N(Wo,z = 2) /dzdWo distribution; see Fig. Qj]), to the 
predictions for the star formation rate at z = 2. Clearly, model 
p redictions a r e in qu alitative agreement with observations. 

IWild et all (120071) have measured the mean star formation 
rate for strong absorbers at intermediate (z ~ 0.6) and high 
(z~ 1.1) redshifts (shown in Fig. [T31 , In particular, all sys- 
tems with Wo > 0.6 were considered in their study. As before, 
we assign a Wo = 0.6A value and note that 68% of those sys- 
tems are expected to lie in the interval 0.6 < Wo < 1.7 A. Our 
mean z = and z = 1 models seem to over predict the data 
by almost an order of magnitude while the median values are 
within a factor ~ 3 of the mea surements. A similar discrep- 
ancy was noted by IWild et aT] (120071) when comparing their 
estima tes f or the SFR to th ose deduced by Prochask a et al.| 
(2005) and iRao et all (120061) . As pointed out by IWild et alJ 
(2007), fiber losses could be important, especially in cases 
where the host extends, or is located beyond, 8kpc from a 
quasar sightline. In fact, our model predicts that much of 
the absorption originates from scales as large as the size of 
the H I disk, and hence beyond the central 8kpc of > 0.01L* 
galaxies. Therefore, the reported values should be consid- 



ered as lower limits on the mean star-formation rates asso- 
ciated with strong systems. Lastly, we note a potential bias 
in probing highly star-forming systems using absorption line 
systems: the Kennicutt-Schmidt law predicts that, on average, 
disks with higher star formation rates would be more optically 
thick, having higher metallicities and, potentially, higher ex- 
tinctions through their disks. It is therefore likely that the high 
SFR end of the distribution of absorbing galaxies will not be 
recorded in surveys where the background quasars are drawn 
from magnitude-limited samples. 

5.7. The Wo — p Relation for Strong Systems 

Our model gives specific predictions concerning the rela- 
tion between Wo and the impact parameter, p, from the center 
of the host galaxy (but not from satellite galaxies that may 
cluster aroun d the true host). In comparing our model pre - 
dict ions to t he Steid el & Sargent! dl992l) . ISteidel et alJ d 19941) . 
and Steidel ( 1995) data set, we assume that these studies iden- 
tify true absorber-host associations. The Steidel et al. sample 
includes the impact parameter between the absorber and the 
host galaxy, the rest equivalent width, and the luminosity of 
the host. To properly compare model predictions to the data, 
we first select a particular absorber-galaxy pair from the Stei- 
del et al. sample. For the chosen system, which is charac- 
terized by a certain galaxy luminosity and impact parameter, 
we calculated the predicted, Wo = Wo(L,p), according to our 
model B prescription, using a Monte Carlo approach. This 
procedure is then repeated for all absorber-galaxy pairs in the 
Steidel et al. sample^ 

The statistical predictions for Wb(p) are compared to the 
observations in Figure [14] in which a typical realization of 
the Monte Carlo simulation is shown. Clearly, the agreement 
is satisfactory in terms of the predicted mean and the scatter 
around the mean. Nevertheless, it is important to emphasize 
that one should be cautious to not over-interpret the outcome 
of this comparison. In particular, to properly account for the 
observations, especially at large impact parameters, one needs 
to account for the possibility of satellite galaxies contributing 
to the absorption. This is currently not included in our model. 
While one does not expect the cross-section of satellite galax- 
ies to fully cover the sky at large impact parameters around the 
host, we note that recent studies show that, in contrast to the 
classic picture, the covering factor for strong Mg II absorption 
in gala xy halos is, on average , smaller than unity on < 50 kpc 
scales (Tripp & Bowen 2005; IChen & Tinker! [2008b . In par- 
ticular, Tripp & Bowen (2005) reported a covering factor of 
order 50% for strong systems. If due to ~ L* hosts then this 
value is qualitatively consistent with the expectations from 
randomly inclined thin disks. A more comprehensive cen- 
sus of absorbers and satellite galaxies on large scales around 
hosts is needed to further test and constrain the model. 

5.8. The Kennicutt-Schmidt Law Revisited 

In §2 we qualitatively showed that strong intervening sys- 
tems obey a version of the Kennicutt-Schmidt law for star 
formation. In particular, by combining observed relations be- 
tween absorbers and star forming galaxies, the power-law de- 

10 We do not show the results for models A and C but note their poor agree- 
ment with the data for the following reason: the Steidel et al. sample suggests 
a decline in the mean Wo for p > 50 kpc. In contrast to that, models A and C 
are characterized by constant wind velocity across galaxy disks. As brighter 
galaxies have larger disks and more intense SFR, these models predict, on 
average, an increase in the mean Wo with increasing impact parameter. 
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FIG. 14. — One realization of our model predictions (see text) for the Wo — p 
relation (black dots) and the data from Steidel et al. (circles). The match 
between model predictions and observations is satisfactory. It should be noted 
that the model does not include the contribution of satellite galaxies to the 
signal. More recent versions of the plot (not shown due to uncertainties on 
the hosts' luminosity) reveal strong absorption around galaxies on even larger 
scales than shown here. 

pendence of the star formation rate on the (absorbing) H I col- 
umn density was deduced. In addition, the correct normaliza- 
tion of the Kennicutt-Schmidt law was obtained when using 
an order of magnitude estimate for the size of the H I disk 
of absorbers' hosts. However, having now obtained a more 
quantitative connection between absorbers and galaxies, the 
proper disk size, as a function of Wo, may be used to arrive at 
a consistent Kennicutt-Schmidt relation for absorbers within 
the framework of the model. More specifically, we need to 
evaluate the expression 



M(Nh)) = 



T[Wo((^V H »] 



7tR eS {L[W {(N H ))}} 2 



(12) 



Aside from Wo((Nn)), all other relations are given by our 
model. In particular, T(Wb) and L(Wq) are drawn from the 
distributions shown in Figure Q whereas R e s(L) comes from 
the fit to the d 2 N(z — 0) /dzdWo distribution, as performed in 
§4 for the different models considered here. 

What Wo((Nh)) relation should we take? The Kennicutt- 
Schmidt law, £*((Nh)), connects the star-formation rate per 
unit surface in the disk to the arithmetic mean of the total hy- 
drogen c olumn density over the disk: (Nh) = (iVm + 2A^h, ) ■ 
For the (lKennicuttlll989[ 119981) samples, it was found that 
(Nm) — 2(A^h 2 ) an d we shall therefore assume that the to- 
tal hydrogen column density, as traced by strong Mg II sys- 
tems, and as average d over a finite patch in the disk, sat- 
isfies (Na) =2{N m ) dFukugita et aljTT998l) . Note that this 
does not necessarily imply that H2 molecules should be seen 
along sight lines to background quasars since galaxy disks are 
highly structured, the H2 gas has a low covering factor, and 
sight lines through molecular gas ma y be obscured by dust 
extinction (Zwaan & Prochaska 2006). 

H I observations of galaxy disks show a very broad rang e 
of column densities over the disk surface dZwaan et al.C 2005), 
and defining the relevant mean, (Nm), is important. As dis- 
cussed in §5.3, the H I distribution for g alaxy disks is con sis- 
tent with that measured for absorbers by Ra o et al.1 (12006ft . In 
particular, for both distributions, the arithmetic mean is much 
larger th an the geometric mean . Fo r strong intervening Mg II 
systems Menard & Chelo uchel (120091) showed that the arith- 
metic mean is about 20 larger than the geometric mean for 
0.5 < W < 0.7 A systems in the dRao et al.ll20 06) sample. In- 



terestingly, in that sample, the measured H I column density 
for strong systems never exceeds 5 x 10 21 cm~ 2 , even though 
the geometric mean over all systems in a given Wo bin does 
increase as, roughly, W 2 (Eq. [TJ. At such a column den- 
sity, the expected reddening for a galactic dust-to-gas ratio is 
E(B — V) ~ 0.3 mag, and the SDSS is expected to miss most 
of the systems resulting in an under-estimation of (Afai), per- 
haps explaining the flat (iVm(Wo)) dependence (Menard et al. 
2008). If galaxy disks are self-similar, as would be expected 
for e.g., a turbulent medium down to the dissipation scale then 
the median and mean would, in fact, be proportional for all 
galaxy luminosities. For these reasons, we shall assume the 
following expression linking the arithmetic mean of the total 
hydrogen column density and WrFT: 



W ((JVh» = 



(Nh) 



1.2 x 10 21 cm- 



0.6 



A. 



(13) 



In doing so we implicitly assume that the iRao et all d2006l) 
data are of relevance to z — systems. This assum ption is sup- 
ported by our findings (§5.3) and the findings of IZwaan et alj 
(2005) concerning the match between the H I cross-section 
of z — galaxy disks and that required to account for high-z 
DLA systems. 

The relevant radius, R e s, to consider in equation [121 is that 
which roughly corresponds to the size of the optical disk 
and over which the Schmidt-Kennicutt law was originally 
evaluated. In our model: R e s(L) ~ riR Mg n(L*)(L/L*y/ 2 ~ 
4(L/L*) 35 kpc, and therefore corresponds, by definition, to 
the core radius discussed in §4. 1 . We ran Monte Carlo simu- 
lations obtaining, as before, a mock absorber-galaxy catalog 
listing (Nh)] for each object, and calculated the mean re- 
lation between the two quantities shown in Figure [TBI 

Comparing our model predictions for the Kennicutt- 
Schmidt law derived for absorbers to the classic Kennicutt- 
Schmidt law for galaxies, we find that the slope of the rela- 
tion matches well the galactic Schmidt-Kennicutt law. The 
obtained normalization has an offset of a factor ~ 1/3 along 
the x-axis or factor ~ 3 along the y-axis compared to the galac- 
tic Kennicutt-Schmidt law (see Fig. fTsT t- Such a discrepancy 
could result from several different reasons: it may be due to 
an under-estimation of the factor between the arithmetic and 
the geometric means of the H I columns in strong systems, 
perhaps as a result of an extinction bias. Alternatively, the 
mismatch could be due to our adopted R e ff(L) being smaller 
by a factor ~ \/3 ~ 1.7 than the true radius. It may also be 
the result of our usage of high-z relations between Wo and Nr\ 
to infer on the prope rties of local star-forming disks [a recent 
theoretical study by Gnedin & Kravtsov (2010) suggests that 
the Schmidt-Kennicutt law may be different at higher-z] . Nev- 
ertheless, given the uncertainties inherent in our analysis and 
the simplifications used, we find the match between the pre- 
dicted Kennicutt-Schmidt law for absorbers and the observed 
one for galaxies to be very good. Furthermore, upon taking, 
for example, the ratio of the geometric to arithmetic mean to 
be 50, we find an excellent match with the Schmidt-Kennicutt 

1 1 We note, for completeness, that we have tried a different set of models in 
which the column density distribution along sight lines piercing t hrough disks 
at cert ain annuli obeys the H I column density distribution of Zw aan et alj 
12005, see their Fig. 3) while the mean column density at any given annuli 
decre ases exponentially with t he impact parameter from the galaxy center 
[e.g., Kalberla & Dedes 1 2008) and references therein]. This model resulted 
in a similar Kennicutt-Schmidt law as the one presented here and will not be 
further discussed due to its more complicated nature. 
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FIG. 15. — Model predictions for the Kennicutt-Schmidt law as deduced 
from model B parameterization while assuming that the disk area over which 
star form ation take s place is given by Tlpg . In constructing the plot we use 
equations 1 13l & [T"2l and plot model B as connected black points. Clearly, the 
predicted Kennicutt-Schmidt relation for absorbers in the qualitative agree- 
ment with the observed Kennicut t-Schmidt law of galaxies [shaded area in- 
cluding uncertainties reported by Kennicutt 1 1998)]. Nevertheless, a small 
correction factor of about factor ~ 2.5 along the *-axis or a factor ~ 0.3 
along the y-axis leads to a better quantitative match with the observations 
(solid black line). All models, A-C, give qualitatively consistent predictions 
(see table [T]for color coding) upon multiplying their N\i values by a factor of 
2.5. Model B's predictions show also the unc ertainty interval in dashed black 
curves. Recent measurements by Biaiel et al. 12008) (thin contour lines) are 
also plotted showing a departure from a simple powerlaw for the Kennicutt- 
Schmidt law for low Nn. Also shown in dotted line is our model predictions 
when the relation between Wq and the geometric mean (instead of the arith- 
metic mean) H I column density is used and assuming (Nu) = (Nri). 

law for galaxies. A qualitatively similar agreement with the 
galactic Schmidt-Kennicutt law is obtained also for models A, 
B, and C (Fig. Q3}. 

To further check the correspondence between absorber and 
galaxi es, it is of c onsid erable interest to extend the Wo range 
of the Ra o et al] (12006b sample and quantify the Kennicutt- 
Schmidt law over a broader range of parameter space. Inter- 
estingly, deviations from the classic Kennicutt-Schmidt law 
at the low Nm end have recently been reported for galaxy 
disks, showing a flattenin g and a subsequen t truncation at still 
lower column densities (Bigiel et al. 2008). Our model pre- 
dictions suggest that most strong, Wo > 0.5A systems arise in 
faint galaxies having relatively low star formation rates. The 
nature of absorption systems is that they probe highly local- 
ized regions in space, so it is of interest to see whether devia- 
tions from a pure powerlaw are consistent between absorption 
and emission derived Schmidt-Kennicutt laws. If such devi- 
ations persist as better statistics are obtained then it may be 
possible to better quantify the correspondence between our 
view of disks in absorption and in emission. This will allow 
to probe faint star-forming disks at high redshifts using ab- 
sorbers' phenomenology in a way which is complementary 
to that employed by emission studies. In particular, these 
approaches suffer from somewhat different systematics and 
their synergy will hopefully help to better quantify star forma- 
tion phenomenon and galaxy evolution physics across cosmic 
time. 

Thus far we have considered a global version of the 
Kennicutt-Schmidt law, as deduced from the phenomenology 
of Mg II systems. Nevertheless, there may also be a local 
version of this law dKennicutt et al.ll2~007l) whereby local out- 
flow properties relate to local star formation processes. At 
present we do not have observational evidence to support this 
hypothesis but note that a relation between the column den- 



sity of O VI ions and the velocity dispersion of their reso- 
na nce 1032A ab s orpt ion transiti on has been recently rep orted 
bv iBowen etafl d2008t see also iHeckman et all (120021) ). In- 
terestingly, the powerlaw slope of this relation is similar to 
the A^hi — Wo relation discussed in this work. It is possible 
that this is a different manifestation of the Kennicutt-Schmidt 
law, highlighting the preponderance of outflows and feedback 
effects from local star forming regions in the Galaxy. 

5.9. IGM Enrichment by Galactic Outflows 

We find no evidence that galactic outflows, as probed by 
strong Mg II systems, extend considerably beyond the H I disk 
of their hosts (i.e., % a ~ 1). This is consistent with a sce- 
nario in which most outflows that give rise to strong Mg II 
systems are launched from compact regions in the disk and 
do not necessarily escape the galaxy and join the IGM. That 
said, our results do not preclude the possibility that galactic 
outflows reach the large scale galactic halo since, on these 
scales, the outflow opacity drops as the gas expands. Map- 
ping Wo > 0.5A systems around nearby galaxies of known 
disk inclinations will shed light on the geometry and origin 
of the gas and help to determine its role in galaxy formation 
and IGM evolution models. In addition, better estimates for 
the spatial distribution of N\\\ > 10 18 cm~ 2 will determine just 
how far above the disk the cool component of galaxy out- 
flows extends. In any case, the connection established here 
between disk outflows and strong intervening systems im- 
plies, by means of number counts, that the phenomenon of 
gas acceleration and ejection from compact star-forming re- 
gions occurs over a broad range of galax y luminosities an d is 
not limited to (bright) starburst galaxies (Heckman 20031 and 
references therein). 

A Chandra imagin g survey of edge-on disks by 
IStrickland et al.l (120041) presented evidence for extended 
X-ray emission at right angles from disks for objects with 
significant star formation rates. These authors found that 
while a galaxy like the Milky Way shows extended X-ray 
emission on a few kpc scales above and below the disk plane, 
(two) fainter galaxies in their sample, with a lower star for- 
mation rate, showed no evidence for such emission. This led 
IStrickland et al.l (120041) to conclude that extended emission 
due to large-scale hot gas bubbles expanding from the disk 
only occurs when T > 0.1 M Q yr _1 . Nevertheless, to explain 
the shape of the d 2 N/dzdWo distribution for Wo > 0.5A, 
our model requires that outflows occur in disks with a star 
formation rate as low as T > 1O~ 2 M yr _1 . Specifically, 
in our model, galaxies with T > 0.1 M yr _1 can account 
for the iNestor et all (120051) data for W > lA. It remains 
to be seen whether the smooth shape of the d 2 N/dzdWo 
distribution (with no apparent break at Wo ~ lA) is due to a 
fortuitous combination of strong- and weak-system popula- 
tions. Alternatively, the absence of hot gas around galaxies 
with low SFRs could indicate that there is little physical 
connection between the cold gas that is traced by Mg II 
systems seen to be outflowing from the disk, and the extended 
X-ray emitting gas extending high above the disk in some 
galaxies. If true, then this may indicate that cold gaseous 
outflows, manifested as strong Mg II systems, are relatively 
compact, although covering most of the disk surface, and 
little can be inferred from their observable properties on 
the ultimate fate of galactic winds. It may also mean that 
while cool outflows are ubiquitous in star-forming systems, 
only in some objects are conditions favorable for launching 
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large-scale X-ray emitting winds dEverett et al.l 120081) . A 
better understanding of galactic winds or fountains may be 
gained by investigating the statistical properties of high- 
ionization line absorbers (Sembach et al. 2003; Grimes et al. 
2009), and throu gh the modeling of energy a n d momentum 
driven outflows (IStrickland & Stevens! 120051 IMurrav et al] 
20051; iDave & Qppenheimerl 120071; iFangano et al] 120071; 
Kobavashi et al]|2007l;lEverett et al.ll2008tlFujita et al] [2009 ) . 

It is worthwhile comparing the outflow velocity of cool 
Mg II outflows to the circular velocity, v c , of galaxies which 
serves as a measure of the escape velocity, v e ~ \Flv c . 
ICourteau et al .1 (120071) have quantified the dependence of the 
circular velocity , v c , of galaxies on their luminosity and find 
that v c oc L° 3 , which is similar to the scaling of the wind ve- 
locity with the SFR. Noting that L* galaxies at z ~ produce 
stars at a rate of ~ 1.7M yr~ , and using equation|2]we find 

-~0.5^, (14) 

Ve 

which implies that for our model A, v w /v e — 0.8 < 1 while 
for our model B, v w /v e ~ 1.8 > 1 within po (see §5.4). We 
emphasize that these values for v w /v e values are maximal in 
the sense that most galaxies will exhibit slower Mg II out- 
flows according to the distribution defined in §3. Therefore, 
if model A is relevant then Mg II outflows probably do not 
escape galaxies and likely form galactic fountains instead. 
Model B is our more favored model since it accounts for the 
observed kinematics of galactic outflows (§5.4), explains the 
dependence of Wo on the impact parameter from the galaxy 
center, p, and is theoretically more plausible (§4.1). In this 
case, judging from the outflow velocity distribution of Figure 
|2j we estimate that about 40% of all galaxies or, alternatively, 
all galaxies for 40% of the time, will eject outflows that could 
escape their potential wells (note that these estimates are sen- 
sitive to uncertain inclination corrections which we ignore in 
the present study). 

We can estimate the mass loss rate that may be associ- 
ated with such outflows assuming they are able to escape the 
galaxy. Specifically, the ratio of the mass loss rate to the SFR 
for a galaxy of luminosity L may be written as 

M{L) „_i 2d 2 N m , _i .... 

^y~Y L WTj^ gn — m p v w ~ ^r^oopc ( 15 > 

where we have taken a simplified disk geometry of radius 
Po = ?7^Mgii and assumed a constant velocity outflow being 
launched perpendicular to it. As the outflow velocity declines 
with p > po, resulting in a wind velocity smaller than the es- 
cape velocity, we neglect the contribution to the mass loss rate 
from radii larger than po. In estimating the mass loss rate, we 
have taken r/, .300pc = r h I (300 pc) = 1 to be th e relevant pres- 
sure scale height of the disk (Ferriere 1998), and a typical 
hydrogen column for the outflowing material of 10 19 cm~ 2 , 
as appropriate for Wo ~ lA systems (lRaoetal]l2006h and 
assuming hydrogen is mostly in H I. We also take a typical 
(one-sided) outflow speed of v w = 70^ w km s (correspond- 
ing to a Wo — 1 .4<j; w A system), with m p being the proton mass. 
We expect equation [15] to have only a very weak dependence 
on galaxy luminosity since Ry^Vw/Y « L°- 7+a3 /L °= const., 
which is consistent with the conclusions of Dalca ntonl (120071) . 
Assuming a 40% duty-cycle for such outflows (see above) 
and relevant £, w > 1 then galaxies expel gas from their cen- 
ters at a rate comparable to the rate at which gas is turned 



into stars. Interestingly, the typical density of the outflow- 
ing cool material used here is Nm/fh ~ 10~ 2 cm~ 3 , somewhat 
lower than the density of the cool ISM phase in the central re- 
gions of our galaxy but in go od agreement with expectations 
on larger, > lOkpc, scales (lFerrierdll998l see her Fig. 2). 
Clearly our estimates are very qualitative as the particular ge- 
ometr y and structure of th e outflowing material is not known 
[c.f. iTremonti et al] d2007l) who interpret similar Mg II out- 
flows in post-starburst galaxies as extending out to lOOkpc 
from galaxy centers in the form of thin shells]. In addition, 
we do not know what fraction of the gas column density is 
bound to the disk and what remainder constitutes the outflow, 
which could lead to our above expression over-estimating the 
mass outflow rate by a factor of a few (assuming solar metal- 
iclity) for a given pressure scale height of the disk. However, 
despite our qualitative estimates and using the phenomenol- 
ogy of strong intervening absorption systems, we arrive at a 
mass loss rate that is of order the star formation rate from 
galaxies, as has been sugge sted by other works based on very 
different sets of arguments (IMurrav et a l. 2005; Bouch l"et al] 
l2007bl:lDaicantonl(2007HDave & Oppenheimeril2007h . 

6. SUMMARY 

We have argued that quasar sight lines that show strong 
Mg II absorption systems (with Wo > 0.5A) in their spec- 
tra, largely probe intervening star forming disks. While a 
considerable fraction of the (H I) column density associated 
with those systems is likely to originate in the gravitationally 
bound star forming disk itself, the rest equivalent width of the 
lines, Wo, is largely set by the outflow kinematics. A version 
of the Kennicutt-Schmidt law is obtained when combining 
the statistical properties of strong intervening systems with 
those of outflowing material from star forming galaxies. In 
this case, the observed cosmic density of intervening systems, 
d 2 N/dzdWo, naturally follows from the luminosity function 
of galaxies and allows us to establish a statistical relation be- 
tween absorbers and their host galaxies. Conversely, knowl- 
edge of the d 2 N/dzdWo distribution and the star-formation 
rate in galaxies may be used to shed light on the properties 
of galactic outflows with preliminary results pointing to an 
agreement between model predictions and kinematic data for 
Mg II outflows. 

The model presented herein gives detailed and testable pre- 
dictions for the association between galaxies and absorbers. 
For example, based on a fit to the number density distribution 
of strong intervening systems, the model suggests that galax- 
ies brighter than > 10~ 2 L* can give rise to strong intervening 
systems, which is in broad agreement with the observations. 
In addition we find that galaxies over a wide range of lumi- 
nosities can give rise to absorbers of a given Wo. More specif- 
ically, Wo ~ lA systems can originate from faint L < 0.1L* 
as well as bright L> L* galaxies, which may help explain the 
lack of clear correlations between galaxy and absorber prop- 
erties in samples with small numbers of quasar-galaxy pairs. 
Detailed predictions for the statistical relations between the 
W) of absorbers, and the star formation rate and luminosity of 
their hosts (and vice-versa) are provided and can be compared 
to observations. The model highlights a potential difficulty in 
the proper identification of the true host galaxies of absorbers 
due to their faintness and the presence of galaxy-galaxy clus- 
tering. This highlights the importance of UV observations 
of nearby intervening systems to properly study the absorber- 
galaxy connection and demonstrates the advantage in probing 



A Kennicutt-Schmidt Law for Intervening Systems 



21 



such the environments at high redshifts using GRBs. Fur- 
ther tests of our model are becoming possible with the recent 
installation of COS on board HST. In particular, it may be 
possible to obtain more reliable estimates for the local cosmic 
density of absorbers, d 2 N(z = 0) /dzdWo, which is crucial for 
appropriately calibrating the model and extrapolating its pre- 
dictions to higher redshifts. Better identification of the true 
hosts of absorption systems and their associated star forma- 
tion rate, particularly for the faintest galaxy hosts, will enable 
us to compare model predictions to observations. Mapping 
the outskirts of galaxy disks and looking for the signatures of 
galactic outflows may lead to a better understanding of galac- 
tic outflows with direct implications for intervening systems' 
phenomenology. 

The simple correspondence advocated here between ab- 
sorbers and star forming disks accounts for several additional 
properties of strong intervening systems. Specifically, by us- 
ing the luminosity-metallicity relation for galaxies and the de- 
duced absorber-disk relation, we can account for the metallic- 
ity of absorbers extrapolated to z = 0. Furthermore, consider- 
ing the evolution in the luminosity function of galaxies with 
redshift, the evolution in the number density of absorbers per 
unit redshift, dN(z)/dz is matched. Interestingly, in its sim- 
plest version, the model suggests that many of the statistical 
properties of star-forming galaxies (such as the specific star 
formation rate, the dependence of galaxy disk radius on lumi- 
nosity, and the outflow velocity dependence on the SFR) have 
not evolved, considerably, since z ~ 2. As model predictions 
are rather sensitive to the scaling laws for galaxies, measure- 
ment of absorber density at higher redshifts may prove to be 
an important probe of galactic physics. In addition, we find 
that model predictions concerning the star formation rate in 
absorbers' hosts are consistent with current observations, yet 
uncertain systematic effects preclude a more comprehensive 
comparison at this stage. In addition, we find that the number 
density of Lya absorbers with Nm > 10 18 cm~ 2 is broadly 
consistent with the results of 21 cm emission studies of lo- 
cal galaxies. Nevertheless, the abundance of absorbers with 
H I columns in excess of 10 213 cm~ 2 seems to be exceed that 
which is inferred locally. At present we do not know whether 
this relates to evolution effects or to the very different beam 
sizes used by each approach. 

Attributing absorbers to gravitationally bound galaxy disks 
and their related outflows, avoids many of the theoretical com- 
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plications associated with models that attempt to identify such 
systems with diffuse halo gas. For example, problems re- 
lated to the formation of gas around galaxies that is cool, 
yet has a high velocity dispersion, are avoided. The need 
to advect metal-rich gas from the centers of galaxies to large 
scales without considerable mixing and disruption due to var- 
ious dynamical instabilities is also alleviated, and the relevant 
physics to consider is just that of the ISM. Cloud stability is- 
sues become irrelevant since the gas is, at least initially, grav- 
itationally bound and a mass reservoir exists to replenish any 
evaporating fraction of it. Our findings point to a different 
physical state (and possibly also an origin) for weak Mg II 
systems whose properties are not accounted for by the model. 

Our results imply that strong intervening systems can be 
used to study the properties of star-forming galaxies at high 
redshift where direct imaging, especially of faint objects, is 
challenging. In particular, our model implies that the evo- 
lution in the density of absorbers across cosmic time closely 
traces the global star formation rate density. Therefore, it may 
be possible to determine the star formation history of our uni- 
verse using absorption line phenomenology. This nicely com- 
plements emission studies that are prone to different system- 
atic effects. A further advantage is that absorption systems 
may be detected to very high redshifts provided enough bright 
background sources, such as quasars and GRBs may be found. 
In addition, metal absorption systems may be used to probe 
the outer parts of galactic disks, estimate their metallicity in 
an independent way, and assist in quantifying the metallicity 
gradient in disks. This has profound implications for galaxy 
formation and evolution, IGM enrichment history, and feed- 
back processes across cosmic time. 
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